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ABSTRACT 


The  resistance  to  flow  caused,  by  artificial  bed  forms 
in  a  flume  was  investigated.  Curves  were  developed  to  relate 
the  roughness  to  the  flow  depth,  and  to  the  dune  shape. 
Roughness  was  expressed  in  terms  of  non-dimensional  forms  of 
Chezy' s  C,  and  the  equivalent  sand  grain  roughness,  kg. 

Tests  were  performed  on  five  different  dune  shapes 
which  had  wavelength  to  height  ratios  of  40,  32,  24,  16,  and  8, 
which  are  of  the  same  order  of  magnitude  as  the  dunes  and  rip¬ 
ples  found  In  nature. 

The  range  of  depth  of  flow  to  dune  height  ratios 
used  in  the  experiments  was  from  4.75  to  2.25,  which  is  also 
comparable  to  the  values  found  under  field  conditions. 

The  effect  of  surface  roughness  was  investigated  by 
coating  the  original  dunes  with  sand,  and  then  repeating  the 
tests.  The  difference  in  resistance  to  flow  due  to  the  surface 
roughness  was  reported. 

The  results  of  these  tests  were  compared  to  results 
of  previous  studies  performed  in  flumes  which  had  mobile  sand 
beds,  and  to  theoretical  formulae  advanced  in  the  literature. 
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CHAPTER  I 


INTRODUCTION 

1. 1  The  Problem 

Many  of  the  rivers  with  which  man  is  concerned 
flow  in  deposits  of  natural  alluvium.  The  slope,,  width, 
depth,  and  trace  of  the  river  are  a  consequence  of  the 
normal  sequence  of  river  flows,  the  characteristics  of  the 
material  in  which  the  channel  is  formed  and  the  type  of 
sediment,  if  any,  transported  by  the  flow.  In  most  rivers 
of  this  description  the  bed  is  in  motion  at  least  during 
peak  flows. 

Depending  on  the  geology  of  the  surrounding  area, 
and  on  that  of  the  regions  upstream,  the  grains  of  the  bed 
material  may  be  classified  as  sand,  gravel,  cobbles,  or 
boulders.  In  a  typical  case  it  might  be  found  that  a  river 
had  a  bed  composed  of  a  mixture  of  different  sizes  of  material 
with  the  greatest  percentage  falling  in  the  sand  size  range. 
This  material  would  be  in  almost  constant  motion,  except  at 
very  low  flows,  with  the  larger  bed  material  moving  only 
during  peak  floods. 

The  sand  grain  particles  move  most  commonly  by 
rolling  or  saltation.  As  a  consequence  of  this  transport, 
during  tranquil  flows,  irregular  ripples  or  dunes  are 
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formed  from  the  bed  material. 

Previous  studies  have  indicated  that  the  resistance 
to  flow  offered  by  the  river  bed  is  a  function  of  both  the 
grain  size  of  the  particles  on  the  surface  of  the  bed.,  and  of 
the  shapes  into  which  the  material  is  moulded.  That  is.,  there 
is  a  frictional  drag  on  the  flow  due  to  the  surface  roughness., 
and  an  additional  resistance  due  to  the  form  drag  of  the  bed 
forms . 

Some  studies  have  roughly  indicated  the  percentage 
of  the  total  resistance  to  the  flow  due  to  grain  roughness,, 
bed  forms  and  other  factors.  However.,  due  to  the  limited 
extent  of  previous  work,,  or  the  preoccupation  of  these  studies 
with  other  facets  of  the  phenomenon.,  little  is  known  about 
the  actual  relation  between  the  form  drag  caused  by  the  dunes., 
and  the  shape  of  the  bed  forms  themselves.  In  fact.,  little 
research  has  been  done  to  estimate  open  channel  roughness ,  or 
resistance  to  flow*  whether  from  natural  or  artificial 
elements . 

It  was  this  lack  of  knowledge  of  open  channel 
roughness  in  general.,  and  a  particular  desire  to  investigate 
the  dependency  of  bed  form  drag  on  dune  shape,  which  prompted 
this  present  investigation. 

1.2  Purpose  and  Scope 

The  purpose  of  this  investigation  was  to  initiate 
a  study  into  the  form  drag  produced  by  bed  forms  in  alluvial 
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channels.  The  present  work  was  conceived  as  a  pilot  project 
to  briefly  investigate  the  practicability  of  such  research, 
to  determine  if  reasonable  results  could  be  obtained  by  the 
methods  used  herein,  and  to  yield  a  rough  idea  as  to  the 
type  of  relation  which  existed  between  the  dune  shape  and 
the  drag  produced. 

In  order  to  facillitate  the  investigation,  the 
research  was  conducted  in  a  laboratory  flume,  using  idealized 
dune  shapes  produced  from  an  artificial  material,  rather 
than  using  a  natural  movable  sand  bed,  or  observing  a  stream 
bed  under  field  conditions.  This  allowed  a  correlation  of 
roughness  results  with  specific  dune  shapes,  and  made 
possible  an  exact  control  of  the  distribution  of  experimental 
data,  rather  than  being  limited  by  the  multitude  of  shapes 
produced  under  normal  conditions. 

In  order  to  make  the  results  of  the  investigation 
comparable  to  field  conditions  the  dune  shapes  chosen  for 
study  were  of  the  same  range  as  those  found  under  natural 
conditions,  in  rivers  or  canals,  and  under  laboratory  conditions 
in  flumes  which  have  mobile  beds.  That  is,  the  wavelength  to 
dune  height  ratios,  L/k,  were  within  the  range  of  values  to 
be  expected  under  natural  conditions. 

In  addition,  the  ratios  of  depth  of  flow  to  dune 
height,  d/k,  used  throughout  the  experiments  were  approximately 
equal  to  those  which  might  be  expected  under  field  conditions, 
or  in  a  mobile  bed  flume. 
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Typical  data  which  show  the  range  of  values  of  L/k 
and  d/k  ratios  under  field  conditions  are  shown  in  Yalin  (l6). 
Similar  data  relating  to  laboratory  flume  experiments  are 
contained  in  Vanoni  and  Hwang  (5).,  or  Guy,  Simons,,  and  Richard¬ 
son  (15). 

The  use  of  artificial  dunes  in  a  laboratory  flume 
simplified  the  problem  by  creating  flow  conditions  that  were 
close  to  those  of  an  ideal  two  dimensional  flow,  and  by 
eliminating  sediment  transport,  and  resistance  to  flow  due 
to  the  meandering  and  channel  irregularities  which  would  be 
found  in  a  natural  situation. 

The  artificial  bed  forms  extended  the  full  width  of 
the  flume,  and  had  straight  crests  perpendicular  to  the  flow. 
This  presented  a  more  regular  pattern  than  a  natural  channel, 
where  discontinuous  and  curved  crests  are  the  rule. 

The  height  of  the  experimental  bed  forms  was  held 
constant  during  the  investigation,  while  the  wavelength  of  the 
recurring  shapes  was  varied.  In  addition,  the  depth  of  the 
uniform  flow  was  varied  throughout  the  testing  procedure  to 
assess  the  change  in  roughness  due  to  the  changes  in  flow 
depth. 

The  initial  inquiry  was  performed  on  smooth  bed  forms. 
Then,  in  an  effort  to  discover  the  order  of  resistance  due 
to  surface  or  grain  roughness,  the  original  dunes  were  coated 
with  sand,  and  the  test  series  repeated,  noting  the  change  in 
the  total  drag. 
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Tests  were  conducted  for  five  different  dune 
shapes,  both  smooth  and  roughened,  giving  a  total  of  ten 
test  series 0  In  each  test  series,  the  depth  of  flow  was 
varied  from  roughly  four  to  ten  times  the  height  of  the 
bed  forms.  In  every  case  the  slope  and  discharge  were 
varied  to  give  the  required  uniform  flow  conditions  at  the 
desired  depth. 

Measurements  of  the  above  parameters,  plus  velocity 
profiles  taken  at  a  number  of  points  along  each  type  of  bed 
form  allowed  calculation  of  Chezy1 s  C,  and  the  Equivalent 
Sand  Grain  Roughness,  ks,  for  each  flow  condition.  Observa¬ 
tion  of  changes  in  these  values  then  permitted  comparisons 
of  the  performance  of  each  dune  shape. 

1 . 3  Review  of  Previous  Research 

A  review  of  the  applicable  literature  on  this  topic 
indicates  that,  as  far  as  can  be  determined,  no  similar 
research  has  previously  been  undertaken  on  rigid  open  channels. 
Numerous  studies  have  been  made  of  the  channel  resistance 
resulting  from  various  types  and  spacings  of  artificial 
roughness  elements,  however,  in  no  case,  did  these  shapes 
approximate  dunes  or  ripples. 

Suryanarayana  (l)  reported  an  investigation  using 
curved  corrugations,  on  a  flume  bed,  at  various  spacings. 

The  results  shown  in  this  work  are  of  the  form: 
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£  =  f  (  los  y  ,  log  y  ] 

/£  \  L  J 


where : 


a 

C 

g 

L 

y 


average  roughness  height 
Chezy's  roughness  coefficient 
acceleration  due  to  gravity 
longitudinal  roughness  spacing 
depth  of  flow 


In  addition,  it  was  reported  that  in  this  study 

Von  Karman's  coefficient,  K ,  was  found  to  vary  from  O.I885 
to  0.19. 

Mirajgaoker  and  Charlu  (2),  and  Chithambaran  and 

Mirajgaoker  (3),  reported  similar  relations  for  roughness 

patterns  consisting  of  natural  stones  placed  at  various  spacings 

In  the  latter  paper  Von  Karman's  coefficient  was  reported  to 
have  remained  fixed  at  O.38. 

Sayre  and  Albertson  (4),  among  others,  have  reported 
Studies  conducted  using  baffles  of  various  sizes,  shapes,  and 
spacings  as  roughness  elements.  They  reported  similar  logarith- 
relations,  and,  m  addition,  introduced  a  parameter  X  to 
describe  spacing.  In  this  study.  Von  Karman's  coefficient  was 
found  to  be  approximately  0.38  and  constant  for  all  spacings. 

The  applicability  of  these  results  to  the  problem 
dealt  with  in  this  present  investigation  is  somewhat  limited, 
however.  There  is  bound  to  be  a  significant  difference  in  the 
values  of  roughness  due  to  isolated  baffles,  rocks,  or  other 
artificial  elements,  and  those  values  due  to  undulations  in 
the  channel  bed  such  as  are  presented  by  dunes  or  ripples. 

These  previous  studies  serve,  however,  to  illustrate  the 
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factors  on  which  total  resistance  depends,  and  to  indicate  the 
rough  trend  of  results  which  may  be  expected. 

A  few  investigations  have  been  done  on  channels 
having  a  mobile  sand  bed,  and  sediment  transport,  in  order  to 
determine  the  amount  of  resistance  due  to  sediment  transport, 
surface  roughness,  and  form  drag  due  to  bed  forms. 


Vanoni  and  Hwang  (5)  indicated  that  the  resistance 
to  flow  caused  by  dune  shapes  was: 

1  =  3.3  log  LR  _  2.3 

/T  TT 


where:  f  =  friction  coefficient  due  to 

the  form  drag 
k  =  average  dune  height 
L  =  average  dune  length 
R  =  bed  hydraulic  radius 

Considerable  data  is  presented  in  this  paper,  both 
from  the  authors'  experiments,  and  from  other  sources.  In 
addition,  data  is  presented  concerning  the  effect  of  the 
surface  roughness  of  the  bed. 


Simons  and  Richardson  (6)  do  not  consider  the  effect 
of  charge  on  resistance,  but  state  that  the  total  roughness, 
C//g,  of  a  channel  may  be  evaluated  by  determining  the  C/7g" 
of  a  plane  bed  of  the  same  material,  and  then  multiplying  by 
a  factor  C*  to  account  for  form  drag.  Their  study  indicates: 

C*  =  f(d,D,S) 

where:  d  =  depth  of  flow 

D  =  particle  diameter  of 
bed  material 
S  =  channel  slope 

In  another  paper  (7)  the  same  authors  do  not  go  so 


. 


. 
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far  as  to  present  a  relation  relating  roughness  to  bed  profiles, 
however,  they  present  a  collection  of  data  from  which  compar¬ 
isons  may  be  made. 

Richardson  and  Simons  (8)  present  equations  relating 
various  bed  profiles  and  channel  resistance: 

for  ripples:  C  =  (7.66  -  Ch3)  log  d  +  0.13  +  11.0 


/g 


v* 


for  dunes: 


where:  R  =  hydraulic  radius  of  flow 
v*=  /gRS  =  shear  velocity 
4RS=  resistance  to  flow  due  to  dunes 
Criteria  are  also  presented  for  predicting  the  type 


of  profile  which  will  form  under  given  conditions.  The  sand 
used  in  part  of  the  research  for  this  paper,  0.93mnij  was 
roughly  similar  to  that  available  for  the  present  investigation. 


Yalin  (9)  presents  a  theoretical  development  of  a 


relation  between  bed  roughness  and  dune  shape,  based  on  the 
head  loss  due  to  the  expansion  and  contraction  of  the  flow 
over  the  bed  undulations.  The  formula  he  presents  is  as  follows: 


where:  kg  =  equivalent  sand  grain 

roughness  of  bed  material 


2 


=  tan 


4>  -  k  f 


2.5  In  11  d 
k 


J 


2L  V 


s 


where:  =  angle  of  repose  of  bed 

material 


and 
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Regarding  the  velocity  distribution  to  be  expected 
in  an  alluvial  channel  with  an  undular  bed,,  Garde  and  Paintal 
(10)  present  data  from  various  sources  which  indicate  that 
Von  Karman's  coefficient,  K  ,  is  not  constant.  It  is  shorn 
that  it  varies  considerably,  depending  on  the  sediment  concen¬ 
tration,  dune  shape,  and  relative  roughness.  If  lots  of  K  against 
these  factors  show  considerable  scatter,  but  indicate  values 
ranging  from  0.15  to  slightly  greater  than  1.0  in  alluvial 
channels . 

Thus,  most  of  the  studies  to  date  have  dealt  with 
artificial  roughness  patterns  of  discrete  elements.  These  are 
not  of  great  value  to  this  work,  except  to  indicate  the  trend 
of  results,  and  important  parameters. 

The  most  valuable  work  in  channel  roughness  has 
concerned  the  resistance  found  in  alluvial  channels.  However, 
only  a  few  attempts  have  been  made  to  directly  relate  rough¬ 
ness  to  form  drag  and  the  shape  of  the  dunes  present,  and  these 
have  been  largely  theoretical. 


CHAPTER  II 


THEORY  OF  OPEN  CHANNEL  FLOW 


2. 1  General 

Due  to  the  complexity  of  the  mechanics  of  fluid  flow, 
especially  that  found  in  open  channel  applications,  many 
simplifying  assumptions  must  be  made  when  proposing  theories 
to  explain  observed  behavior.  A  standard  assumption,  commonly 
used,  is  that  steady,  two-dimensional  flaw  exists. 

However,  numerous  studies,  such  as  those  of  Keleugan 
(ll),  Tracy  and  Lester  (12),  or  Gruff  (13),  have  demonstrated 
that,  due  to  sidewall  effects,  flow  in  an  open  channel  may  be 
considered  uo  approach  the  ideal  two-dimensional  conditions 

only  in  a  relatively  narrow  region  about  the  channel  centre 
line . 

Cruff  demonstrates  that  within  the  width  to  depth 
ratios  to  be  used  in  this  investigation,  0.111  to  0.246, 
the  shear  stress  distribution  on  the  channel  bed  approaches 
a  uniform  condition  only  in  a  region  about  the  centre  line  of 
the  channel,  approximately  ten  percent  to  twenty  percent  of 
the  channel  width. 

Thus,  it  should  be  carefully  noted  that  implicit 
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in  the  discussion  that  follows  is  the  assumption  that  measure¬ 
ments  are  to  he  taken  within  a  suitable  region  in  which  the 
flow  conditions  can  be  considered  to  approximate  those  of  a 
two-dimensional  flow. 

2.2  Chezy* s  Relation 

The  uniform  flow  relation  developed  by  French  engi¬ 
neer  Antoine  Chezy  about  1769  was  probably  the  first  such 
formula.  It  is  still  widely  used,  the  empirical  constant  being 
estimated  through  experience.  No  suitable  means  has  been 
found  for  calculating  the  constant  from  roughness  measurements 
or  other  factors.  The  observations  of  this  study  may, however , 
present  a  somewhat  clearer  idea  of  the  manner  in  which  certain 
factors  affect  the  final  result. 

Chezy  based  his  mathematical  derivation  on  two 
assumptions : 

1)  That  the  force  resisting  the  flow.,  per  unit  s,rea 
of  stream  bed.,  was  proportional  to  the  square 

of  the  mean  velocity,  and 

2)  That  in  uniform  flow  the  resistance  of  the  bed 
is  exactly  equal  to  the  component  of  the  gravity 
force  which  causes  flow  down  the  slope  of  the 
channel . 

The  Chezy  relation  is  usually  expressed  as  follows: 

v  =  cJrs 

where:  V  =  mean  velocity  of  flow 

C  =  Chezy*  s  constant 
R  =  hydraulic  radius  of  flow 
S  =  slope  of  the  energy  line 


. 

•  | 
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The  formula  may  be  derived  in  the  following  manner: 
(see  fig.  2.1) 


Fig.  2.1:  Development  of  Chezy1  s  Equation 


The  total  weight  ¥  of  the  flow  between  the  two 
sections  shown  is: 


¥  =  r  L  A 


where:  y  =  unit  weight  of  flow 


L  =  length  of  section 
A  =  cross  section  area  of  flow 


The  component  of  the  weight  which  is  causing  flow 
down  the  slope  is: 

¥  sin  S  =  AS 

where:  S  =  slope  of  energy  line  or 

channel 

If  we  assume  a  unit  resisting  force  7'Q  on  the  bed 
of  the  channel  we  find  the  total  resisting  force  to  be: 

X  L  P 


where:  P  =  wetted  perimeter  of  flow 


. 
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Now,  using  Chezy* s  assumption  that  accelerating  and 
resisting  forces  are  equal,  and  assuming,  in  addition,  that 
S  is  small,  we  have: 

y  L  A  sin  S  =  T  L  P 

o 

or:  /AS  =  1  P 

o 

r0  =  A  .S  =  RS  (2.1) 

r  p 

Now,  from  Chezy* s  other  assumption,  we  may  write: 

X  =  K  V2 

where:  K  =  constant  of  proportionality 

V  =  mean  velocity  of  flow 

Therefore:  K  V^  =  RS 

v  =  rv .  f  rs 

/  k 

Replacing  the  constants  y'  and  K  by  a  single  constant 
named  "Chezy1 s  C"  we  write: 

V  =  c/RS  (2.2) 

Thus  the  Chezy  relation  relates  the  mean  velocity 
of  flow,  the  channel  geometry,  and  the  slope.  The  constant 
is  dependent  upon  fluid  properties,  and  the  shape,  distribution, 
and  size  oi  the  roughness  elements  on  the  channel  boundary. 

2.3  Velocity  Distribution 

It  is  possible  to  develop  an  expression  for  the 
velocity  distribution  in  a  steady,  uniform,  two-dimensional 
flow  having  a  free  surface.  Then,  using  this  type  of  expres¬ 
sion,  a  relation  may  be  obtained  between  the  bed  roughness 


. 
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and  flow  parameters. 

From  the  previous  development  of  the  Chezy  equation 
we  obtained  the  following  relation: 

T0  =  rs  r  (2.1) 

Generalizing  this  expression  for  any  depth  within 
a  flow.,  and  for  the  case  of  a  wide  channel.,  we  could  write: 

T  =  y  s  a 

which  would  generate  a  shear  stress  distribution  as  shown  in 
fig.  2.2: 


Fig.  2.2:  Shear  Stress  Distribution  in 

Open  Channel 


At  any  level  y  above  the  bed  the  total  shear  stress 
T  is  composed  of  a  viscous  component,  and  a  turbulent 

component,  Tf 

The  viscous  shear  stress  is: 


r 


i 


=  dv 


dy 


where  dv/dy  is  the  velocity  gradient:  at  the  level  y  under 
consideration. 
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The  turbulent  shear  stress  is: 

0"  t  =;p/>  u'  v' 

where  u' v;  is  the  long  term  average  product  of  the  lateral  (u) 
and  transverse  (v)  components  of  the  turbulent  velocity  fluc¬ 
tuations.  Thus,  for  a  turbulent  flow: 

■y  =  h  +  rt 

However,  as  can  be  seen  from  fig.  2.2,  the  contribu¬ 
tion  of  the  viscous  shear  forces  is  significant  only  within  a 
small  region  relatively  close  to  the  bed. 

Therefore,  at  any  appreciable  distance  above  the  bed 
we  could  write: 

r  *  rt 

As  an  analogy  to  laminar  flow  we  could  write: 


where :  x't  =  turbulent  or  eddy  viscosity 
and  where y  ^  varies  as  some  function  of  y.  This  variation 
of  ye  £  with  y  is  not  known,  but  we  may  say  as  we  approach  the 
boundary: 

lim  u.  =  (const)/  v*  y  (2.5) 

(y-0K  ^ 

where:  v*  =  shear  velocity 

Now,  as  was  previously  mentioned,  the  assumption 
'Y  ^  is  only  valid  above  some  certain  depth  in  the  flow, 

that  is,  say  when  y  ^  <f  ,  where  /  is  the  lower  limit. 

In  addition,  expressions  (2.4)  and  (2.5)  are  only 
valid  in  a  space  completely  occupied  by  fluid,  and  can  not  be 


’ 
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valid  in  the  region  of  roughness,,  0<y<k  where  space  is 

s 

occupied  partly  by  fluid  and  partly  by  boundary  material. 

Thus  ks  also  constitutes  a  lower  limit  of  validity. 

Therefore we  can  say  that: 

y*t  =  ( const )^3  v*  y  (2.6) 

is  valid  if: 


d  ^  y  >/  1 

where:  1  =  larger  of  £  or  k  . 

s 

y  =  elevation  in  question 
Now  when  y -*■(),  then  y  -*■'!' So,  considering  that 

2 

%  =  v*. -  then  by  combining  equations  (2.6)  and  (2.4)  we 

obtain: 

v*  =  ( const ) . y . dv  (2.7) 

,  .  dy 

Integrating  (2.7)  in  the  region  1  to  y,  we  obtain: 

v  -  Xa  =  (const)ln  y 
V*  v*  1 

where:  v^  is  the  value  of  v  at  y  =  1. 

Now,  if  /  >  k  ,  i.e.  if  1=  £  ,  we  have : 

s 

v  _  (const) In  y  +  vs 

V*  /  v*  (2.8) 

and  if  kg  >  /  j  i.e.  if  1  =  kg  we  have: 

v  =  (const) In  y  +  v*f  (2.9) 

v*  kg  v* 

Both  expressions  (2.8)  and  (2.9)  can  be  simplified 

into  the  form: 

(const)  In  y  +  B 
k 


v 

V* 


(2.10) 
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where:  B  =  (const) In  k  +  (const)  if  k/c -*  0  (2. 10)' 

B  =  ( const  )^  if  k  /s—  oo  ( 2 . 10 ) y/ 

Expressions  (2.10)  and  (2.10)  may  be  shown  in  terms 
of  v*  and  v  instead  of  £  as  follows: 

B  =  (const) In  v*kg  +  ( const )  if  v*kg  0 


v 

B  =  (const), 

•“■S 


if  v*.k_ 
*  s 

zo 


CD 


Fig.  2.3:  Experimental  Values  of  B 


The  form  of  experimental  data  is  shorn 
Using  this  data,,  the  following  can  be  deduced: 

B  =  2.5  In  v*k  +5.5  if 

O 

V 


in  fig.  2.3. 


v*ks  <  5 

7? 


B  =  8.5 

N 

and  equation  (2.10)  now  becomes: 

v  =  2.5  In  y  +  B 

v*  k 

s 


lf  Pfs  >  70 

z> 


(2.11) 
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The  case  when  (v*^s/zA  )<  5  corresponds  to  hydraulicaly 
smooth  flow  where  The  case  where  (v*k  A>  )  >  70  corre- 

sponds  to  fully  developed  turbulent  flow.  In  between  these 

cases,  the  flow  Is  transitional  and  values  of  B  must  be 

estimated  from  the  experimental  graph. 

It  is  useful  to  note  that  the  expression  (2.11)  may 
be  simplified  to: 

v  =  2.5  In  A-l  y  (2.12) 

v*  k 

s 

.  „  B/2.5 

where:  A-^  =  e 

For  purposes  of  this  investigation,  the  above  expres¬ 
sion  was  solved  for  kg  yielding: 

X 

ks  =  y-e  (2.13) 

where:  x  =/  B  v  \ 

(23  27333 

Thus,  using  this  relation,  the  equivalent  sand  grain 
roughness  of  any  given  boundary  configuration  may  be  found 
if  the  viscosity,  slope,  depth  of  flow,  and  flow  velocity  at 
a  given  elevation  above  thd  bed  are  known. 

In  a  similar  manner  it  is  possible  to  derive  an 
equation  relating  the  equivalent  sand  grain  roughness  to  the 
mean  velocity  of  the  flow: 

jf  =  2.5  In  11.0  d  (2.14) 

v*  ks 

Equation  (2.11)  may  also  be  written  in  the  following 


manner: 


' 


. 
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v  _  5-75  log  y  +  B 

v*  ks 

With  regard  to  the  assumption  of  two-dimensional 
flow  in  the  central  portion  of  a  channel,  Tracy  and  Lester  (12) 
have  shown  that  some  variation  may  be  found  in  the  factor  5.75, 
which  is  really  2.3 /k,  where  K  is  Yon  Karman's  coefficient, 
depending  on  the  width  to  depth  ratio  of  the  channel.  The 
possible  variation  of  this  factor  in  the  present  study  will 
be  discussed  in  a  later  section. 

2.4  Dimensional  Considerations 

Since  expressions  relating  the  flow  properties  have  been 
developed,  it  is  appropriate  to  examine  the  relationships 
between  the  parameters  involved,  in  order  to  devise  methods 
of  experimental  investigation.  By  varying  various  parameters 
while  holding  others  constant,  it  is  possible  to  devise  appro¬ 
priate  experimental  procedures  which  will  yield  information 
concerning  the  physical  laws  underlying  the  behavior  of  the 
system  under  study. 

The  mean  velocity  observed  in  a  steady  uniform  flow 
can  be  said  to  depend  on  both  the  properties  of  the  fluid,  and 
those  of  the  channel.  In  general,  we  could  write: 

v  =  fi( ksJ  gS  ) 

where :  =  unit  mass  of  fluid 

y*  =  dynamic  viscosity  of  fluid 
d  =  flow  depth 

k  =  equivalent  sand  grain  roughness 
s  of  channel 

g  =  acceleration  of  gravity 
S  =  slope  of  channel 


. 
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Note  that  in  the  above  expression  there  is  no  term 
describing  the  channel  width  or  shape.  This.,  therefore.,  assumes 
that  the  flow  considered  is  in  a  wide  channel,  where  it  may  be 
considered  two-dimensional,  and  unaffected  by  any  sidewall 
effects. 

Here  we  have  used  one  term, k  ,  to  describe  the  ulti- 

s 

mate  channel  roughness.  However,  this  term  may  be  broken  down 
as  follows: 

ks  =  f2(L,k>  jzUx)  (2.15) 

where:  L  =  length  of  bed  forms 

k  =  height  of  bed  forms 
X  =  some  measure  of  surface 
roughness 

-  angle  of  repose  of  bed 
material,  or  slope  of 
dune  face 

Again,  we  have  used  one  term,  X,  to  cover  all  the 
effects  of  surface  roughness.  This  term  could,  in  turn,  be 
broken  down,  and  replaced  by  others  describing  the  roughness 
spacing,  height,  shapes,  size  distribution,  and  other  factors. 
However,  since  in  this  investigation  we  will  not  consider  the 
manner  in  which  roughness  depends  upon  these  characteristics, 
we  may  drop  the  parameter  X  from  our  analysis. 

Similarly,  the  parameter  $  is  constant  throughout 
our  tests,  and  need  not  be  considered.  Therefore,  we  may  write: 

kg  =  f3(L,  k) 

Combining  this  with  our  previous  expression,  we 

obtain: 

V  =  j  d,  L,  k,  gS  ) 
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Rearranging  these  terms  to  form  non-dimensional 
parameters  yields: 


V 

r  gds 


fc  /  ^  L’  VgdS  .  d 
(  k  k  ^7> 


V 

v* 


_  fr 


v*d 


■g  f*>  b 
D  (  k  k  i; 

If  it  is  assumed  that  in  these  investigations  fully 
developed  turbulent  flow  will  exist,  then  the  term  v*d/2; 
may  be  ignored.  Thus,  we  have: 


V 

v 


'6 


(a3  l 
l  k  k 


now,  recall  that: 


V 

v 


* 


Therefore,  the  expression  becomes: 

f, 


<=(1-0 

Again,  equation  (2.15)  yields: 

f 1(L,  k,  0  ,  X) 


C 

/s 


k 


(2.16) 


or,  we  could  write: 


k. 


Therefore : 


n.k 


where:  n  =  f  (k,  L,  0>  ,  X) 


kg  =  f2(k,  L,  ft  ,  X)k 


N  Again,  since  they  are  constant  in  this  study,  the 
parameters  <f>  and  X  may  be  neglected.  Also,  as  before,  it 
will  be  noted  that  the  presence  of  a  wide  channel  is  assumed 


. 
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since  no  channel  geometry  terms  are  included. 

The  expression  therefore  reduces  to: 
ks  =  f3(k,  L)k 

Therefore: 


k 

0s  = 

or:  k 

s 

Rearranging  these 


f4(k,  ■  L)k 
d 

f5(r  k,  d) 

terms  gives: 


(2.17) 


Examination  of  relations  (2.15)  and  (2.17)  suggests 
forms  of  results  which  would  most  appropriately  demonstrate 
the  laws  governing  the  resistance  of  bed  forms  to  uniform  flow. 


\ 
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CHAPTER  III 


THE  EXPERIMENTAL  INVESTIGATION 

3. 1  The  General  Approach 

The  testing  procedure  consisted  of  placing  artificial 
bed  forms  in  a  flume,  and  establishing  uniform  flow  over  them. 
Velocity  profiles  were  taken  at  the  centre  line  of  the  flume 
at  a  number  of  sections.  In  addition,  the  depth  of  flow, 
discharge,  and  flume  slope  were  recorded. 

Tests  were  performed  for  bed  forms  having  wavelength 
to  height  ratios  of  40,  32,  24,  l6,  and  8.  For  each  shape 
tests  were  run  for  depths  of  flow  ranging  from  4.5  to  9.5 
times  the  amplitude  of  the  bed  forms. 

An  initial  test  series  was  run  using  smooth  concrete 
bed  forms.  These  bed  forms  were  then  roughened  by  adding  a 
coat  of  sand,  and  the  test  series  was  repeated. 

By  using  the  data  observed  during  these  tests  it  was 
possible  to  calculate  the  roughness  parameters  as  outlined 
previously,  and  to  determine  the  manner  in  which  this  char¬ 
acteristic  roughness  varied  with  the  bed  form  shape  and  the 
depth  of  flow.  Details  of  the  equipment  used,  and  the  procedure 
followed  are  given  in  the  following  sections. 
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3. 2  The  Test  Flume 

The  flume  used  for  this  investigation  was  of  the 
recirculating  type,,  having  its  own  reservoir  and  pump.  In 
addition,,  provision  was  made  for  the  addition  or  release  of 
water  from  the  system  to  a  larger  sump  within  the  building. 
Thus.,  although  each  day' s  operation  was  normally  performed 
using  only  the  water  in  the  flume  system.,  it  was  possible  at 
any  time  to  discharge  this  water  into  the  larger  sump  where 
it  could  be  filtered,  or  to  add  fresh  water  from  the  larger 
system  to  make  up  for  evaporation  or  leakage  losses  in  the 
flume . 

The  flume  was  of  aluminum  frame  construction,  having 
plexiglass  walls,  a  teak  tailgate,  and  a  steel  headbox,  reser¬ 
voir,  and  piping.  It  was  pivoted  at  the  headbox  outlet,  and 
had  a  vertical  screw  adjustment  near  the  tailgate  to  allow 
simple  adjustment  of  the  slope.  The  slope  was  read  directly 
from  a  graduated  scale  at  this  point  on  the  side  of  the  flume. 

The  inside  dimensions  of  the  flume  were  as  follows: 
width  19.625  inches,  depth  29.0  inches,  length  16.25  feet. 

A  layout  diagram  of  the  flume  is  shown  in  fig.  3.1. 

The  downstream  gate  was  formed  of  two  sets  of 
vertical  fingers  which  could  be  raised  independently,  allowing 
reasonably  easy  control  over  the  depth  of  flow,  although  some 
warping  of  the  wood  caused  occasional  sticking  of  the  gate. 

A  one  piece  aluminum  undershot  gate  was  provided 
at  the  headbox  outlet,  but  this  was  not  used  during  the  tests. 
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FIG.  3-1  EXPERIMENTAL  EQUIPMENT  -  GENERAL  LAYOUT 
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Plate  3.1  :  General  Layout  of  Laboratory  Equipment 


Plate  3.2 


Mini-Flowmeter  and  Counter  Unit 
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Plate  3.3  :  Bed  Forms  and  Flowmeter  in  Flume 
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Control  of  the  discharge  was  accomplished  by  a  valve 
downstream  of  the  constant  speed  centrifugal  pump. 

Expanded  metal  mesh  was  installed  at  the  entrance 
to  the  headbox  to  prevent  concentration  of  the  flow.  In 
addition,  during  the  test  series  conducted  with  the  smooth 
concrete  bed  forms,  about  the  first  four  feet  of  the  flume  bed 
were  covered  with  the  same  material  to  aid  in  the  quick 
development  of  fully  developed  turbulent  flow,  and  the  estab¬ 
lishment  of  a  logarithmic  velocity  profile. 

To  effect  small  changes  in  the  tailgate  setting  it 
was  occasionally  found  to  be  better  to  place  similar  sheets 
vertically  in  front  of  the  tailgate  than  to  tamper  with  the 
gate  setting.  The  slight  reduction  in  flow  area  thus  achieved 
served  for  minor  adjustments  in  the  flow  depth. 

It  was  found  that  after  standing  idle  for  a  number  of 
days  considerable  rust  sediment  collected  in  the  system  pipes. 
This  residue  was  deposited  on  the  plexiglass  sides  of  the 
flume,  and  on  the  velocity  measuring  probe,  and  contributed 
an  unpleasant  odor.  In  such  a  case  the  flume  system  was 
repeatedly  flushed  with  clear  water  from  the  main  building 
sump,  and  the  waste  water  filtered.  The  build-up  of  such 
deposits  occasionally  necessitated  the  cleaning  of  the  velocity 
probe  and  the  flume  sides. 

The  general  arrangement  of  the  flume  is  shown  in 

Plate  3.1. 

The  discharge  was  measured  by  means  of  a  magnetic 
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flow  meter  installed  just  downstream  from  the  pump.  This  device 
produced  an  electrical  potential  proportional  to  the  discharge,, 
and  was  attatched  to  a  sensitive  voltmeter.  For  each  test 
run  the  reading  of  the  voltmeter  was  recorded,  and  the  discharge 
was  read  from  a  calibration  curve; 

Depth  of  flow  was  measured  by  means  of  two  trans¬ 
parent  scales  affixed  to  the  sides  of  the  flume.  The  water 
surface  in  the  flume  could  be  observed  through  the  scales  and 
ohe  flume  sides,  and  the  depth  read  directly.  During  the 
course  if  the  investigation  all  depths  were  measured  relative 
to  the  minimum  point  of  the  bed  form  profile.  However,  in 
order  to  allow  a  more  realistic  approximation  to  natural  river 
conditions  the  datum  for  all  elevation  and  depth  measurements 
was  later  moved  to  a  point  midway  between  the  high  and  low 
points  of  the  dunes.  That  is,  the  average  elevation  of  the 
bed  was  adopted  as  the  elevation  datum.  The  tables  of  obser¬ 
vations  in  Appendix  A,  and  all  plots  of  results  are  based  on 
this  datum. 

3*  3  The  Velocity  Meter 

The  velocity  profiles  were  measured  by  a  Gloster 
Mini-Flowmeter,  and  a  Dekatron  Counter  Unit.  The  flowmeter, 
or  probe,  was  1  cm.  in  diameter,  and  hence  it  was  possible  to 
obtain  velocities  relatively  accurately  at  a  number  of  eleva¬ 
tions  above  the  bed.  Readings  were  taken  at  half-inch  inter¬ 
vals,  generally  from  an  elevation  of  three  quarters  of  an  inch 
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above  datum ,  to  within  one  quarter  of  an  inch  of  the  surface. 

In  a  number  of  cases  it  was  found  that  when  the 
probe  was  very  close  to  the  water  surface,  say  within  one 
inch,,  the  readings  obtained  were  somewhat  less  than  expected. 
This  was  due  to  the  fact  that . the • flowmeter  was  designed  to 
operate  under  a  certain  depth  of  water  to  ensure  proper  func¬ 
tioning.  When  the  depth  over  the  probe  was  less  than  this 
specified  depth,  the  propeller  occasionally  did  not  register 

the  correct  number  of  revolutions  for  the  incident  velocity. 

/ 

;  This  defect  was  overcome  by  plotting  each  velocity 

profile  taken.  In  the  case  of  those  profiles  which  showed 
abnormal  curvature  near  the  water  surface  the  top  portion 
of  the  profile  was  disregarded,  and  observed  values  within 
that  region  were  not  used  in  subsequent  calculations.  Those 
values  so  neglected  are  indicated  by  an  asterisk  in  the  data 
tables  in  Appendix  A. 

It  was  found  that  unless  the  water  in  the 'flume 
system  was  exceptionally  clear  a  coating  of  rust  residue  was 
deposited  on  the  propeller  and  guard  of  the  probe.  In  order 
to  avoid  any  retarding  of  the  meter  from  this  effect,  it  was 
necessary  to  wash  the  probe  in  a  weak  solution  of  Hydrochloric 
acid,  and  to  rinse  it  in  distilled  water.  This  was  generally 
done  at  the  end  of  each  day’s  test  run. 

In  addition,  slight  trouble  was  encountered  due  to 
the  tendency  of  the  probe  to  collect  bits  of  hair  or  lint 
which  were  present  in  the  water.  Any  build-up  of  such  material 
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tended  to  foul  the  propeller,  and  necessitated  a  tedious 
cleaning.  However,  in  general,  the  ease  of  operation  of  the 
flowmeter  unit  outweighed  any  disadvantages,  and  this  method 
proved  most  satisfactory  for  velocity  measurements. 

Some  electrical  trouble  was  encountered  with  the 
Dekatron  Counter  unit  but,  in  general,  it  provided  a  satis¬ 
factory  method  for  controlling  the  flowmeter. 

It  was  noticed,  however,  that  the  length  of  coaxial 
cable  used  to  connect  the  probe  to  the  counter  unit  had  an 
effect  on  the  number  of  revolutions  registered.  In  order  to 
minimize  this  resistance  the  shortest  possible  length  of 
cable  was  used  throughout  virtually  all  tests. 

The  flowmeter  and  counter  unit  used  are  shown  in 
Plate  3.2,  and  may  be  seen  In  operation  in  Plates  3.1,  and  3.3. 

A  recallibration  was  performed  on  the  mini-flow¬ 
meter  after  the  laboratory  tests  were  completed.  The  veloci¬ 
ties  calculated  from  the  original  manufacturer's  calibration 
curve  were  found  to  vary  from  six  percent  low  at  velocities 
of  one  half  foot  per  second,  to  three  percent  low  at  velocities 
of  one  and  one  quarter  feet  per  second. 

3.^  The  Bed  Forms 

The  bed  forms,  or  dunes,  were  cast  in  concrete  as 
shown  in  fig.  3-2.  The  height  of  the  individual  dunes  was 
held  constant  at  one  half  an  inch,  and  the  wavelength  was 
varied  to  produce  the  desired  shapes.  The  dunes  were  cast 
either  singly,  or  in  multiple  units,  depending  on  their  size. 
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FIG.  3-2  BED  FORM  PROFILES  AND  CONSTRUCTION 
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so  as  to  provide  easily  handled  units. 

The  shape  of  the  dunes  was  somewhat  idealized  to 
facilitate  fabrication.  All  surfaces  on  the  dune  were  plane 
surfaces,  not  curved  as  would  be  found  in  nature.  However,, 
it  was  felt  that  for  the  purpose  of  the  investigation,  and 
within  the  accuracy  of  the  experiment,  little  actual  error 
would  result  from  this  simplification. 

The  lee  slope  of  the  dunes  was  chosen  as  38.6.°  This 
angle  was  chosen  since  it  was  easy  to  manufacture  when  making 
thp  forms  for  the  concrete  dunes  (rise  over  run  equal  to  4/5) 
and  since  it  was  close  to  the  natural  angle  of  repose  of  sub¬ 
merged  sand.  This  angle  was  used  in  all  dune  shapes. 

The  width  of  the  dune  modules  was  chosen  as  19.5 
inches.  This  allowed  a  reasonably  tight  fit  in  the  flume, 
having  only  l/l6  of  an  inch  on  each  side,  but  allowed  easy 
installation  and  removal  of  the  bed  form  blocks. 

The  dune  modules  were  simply  placed  on  the  floor  of 
the  flume,  and  held  in  place  by  their  own  weight.  This  proved 
entirely  satisfactory,  and  no  movement  or  shifting  of  the 
blocks  was  observed  at  any  time. 

The  actual  length  of  the  portion  of  the  flume  bed 
covered  by  the  bed  form  blocks  varied  slightly  from  one  test 
series  to  the  next,  depending  on  the  integral  number  of  blocks 
which  most  closely  approximated  the  length  of  the  flume. 
However,  in  no  case  was  the  covered  portion  less  than  the 
16 1  -3,r  distance  from  entrance  gate  to  tailgate.  Any  excess 
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block  length  was  allowed  to  project  into  the  headbox. 


After  the  first  test  series  using  the  smooth  concret 
dunes,  each  set  of  bed  forms  was  coated  with  sand,  and  the 
test  runs  were  repeated.  The  sand  used  was  white  Ottawa  sand 
having  an  average  diameter  of  0.0282  inches,  or  0.716  mm.  A 
description  of  a  grain  size  analysis  performed  on  this  mater¬ 
ial  is  contained  in  Appendix  A. 

The  procedure  followed  in  coating  the  dunes  with 
sand  was  as  follows:  The  dunes  were  first  given  a  light  coat 
of  varnish,  and  allowed  to  dry.  Each  module  was  then  heavily 
coated  with  varnish  and  was  buried  in  a  box  filled  with  the 
Ottawa  sand.  The  varnish  was  allowed  to  set  for  twenty  min¬ 
utes,  then  the  bedform  block  was  carefully  unearthed,  and 
allowed  to  dry  overnight.  Ail  excess  sand  was  then  brushed 
off,  and  any  bare  spots  were  patched  with  a  light  coat  of 
varnish  and  sprinkled  sand.  At  the  same  time,  any  lumps  on 
the  surface  were  smoothed  out  as  well  as  possible.  When  all 
the  blocks  were  dry  and  evenly  coated  that  set  was  deemed 
ready  for  its  repeat  test  series. 

The  final  appearance  of  the  blocks  after  coating 
was  of  a  rough,  but  generally  uniform, layer  of  sand  usually 
two  to  three  grains  thick.  The  contours  of  the  original  bed 
forms  were  closely  followed,  and  there  was  no  noticeable 
tendency  for  the  sand  to  be  thicker  in  the  low  points  of  the 
bed  profile.  That  the  varnish  was  reasonably  waterproof  is 
shown  by  the  fact  that  very  little  sand  was  observed  to  wash 
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off  the  dunes.  The  sprinkling  that  was  found  in  the  flume 
reservoir  at  the  end  of  the  investigation  was  most  likely 
loose  sand  that  had  not  been  brushed  off  the  blocks  before 
testing,  or  sand  knocked  loose  during  installation  or  removal 
of  the  dune  modules  from  the  flume. 

3*5  The  Experimental  Procedure 

The  test  procedure  for  all  runs  was  as  follows: 

The  chosen  bed  form  blocks  were  set  into  the  flume,  and  the 
depth  scales  zeroed  on  the  low  points  of  the  bed  form  profile. 

A  desired  depth  of  flow  was  then  chosen,  and  the  discharge, 
tailgate  setting,  and  the  slope  of  the  flume  were  adjusted 
until  uniform  flow  was  achieved  at  the  desired  depth.  The 
discharge  and  slope  for  that  run  were  then  noted.  Using  the 
mini -flowmeter  probe  a  velocity  profile  was  .then  taken  at  the 
centre  line  of  the  flume  at  three  different  sections  along  a 
dune.  (See  fig.  3.3.) 

m  order  to  obtain  velocity  profiles  observations 
were  taken  at  increments  01  one  half  inch  above  the  observation 
datum.  At  each  point  an  average  was  taken  from  five  ten-second 
counts  of  the  probe  revolutions.  Thus,  in  effect,  each  value 
was  the  mean  over  fifty  seconds  of  the  velocity  at  that  point. 

This  procedure  was  repeated  for  depth  increments 
of  one  quarter  inch,  from  two  and  one  half  inches  to  five 
inches  above  the  observation  datum. 

Bjy  using  the  ocservations  recorded  the  required  results 
were  then  determined  by  calculation. 
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It  should  be  noted,  again^  that  a  change  was  made 
at  this  point  between  the  datum  used  for  experimental  obser¬ 
vations.,  and  the  datum  to  which  all  calculations  and  graphs 
are  related. 

For  ease  of  operation  the  elevation  datum  during  all 
test  runs  was  taken  as  the  low  point  of  the  dune  shape  profile. 
However,  during  the  calculations  one  quarter  inch  was  subtracted 
from  all  depths  of  flow  and  elevations  above  the  bed,  so  as 
to  move  the  datum  to  a  level  equal  to  the  average  height  of 
the  bed  profile.  This  is  shown  in  fig.  3.3. 

It  was  felt  that  this  new  datum  was  a  more  natural 
one  If  the  experiments  were  to  be  compared  to  natural  field 
conditions.  In  addition,  if  the  tests  were  to  be  repeated 
using  a  mobile  sand  bed,  then  this  new  calculation  datum  would 
equal  the  original  elevation  of  the  plane  bed  at  the  start  of 
the  test,  after  which  dunes  would  form,  the  amount  of  material 
above  this  datum  being  equal  to  the  amount  required  to  fill 
the  ''valleys’1  below  the  datum. 

In  the  calculation  of  the  results,  first  all  velocity 
profiles  were  plotted,  and  low  values  of  velocity  near  the 
surface  were  eliminated,  as  explained  in  section  3*3*  Then, 
using  equation  (2.13)^  individual  values  of  ks  were  calculated 
from  each  velocity  reading  taken  in  the  velocity  profiles. 

In  performing  this  calculation  it  was  first  assumed 
that  (v*k s/jj  )  would  be  greater  than  JO,  yielding  a  B  value 
of  8.5.  Using  this  value  of  B  the  calculation  was  then  per- 
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formed,  and  the  value  of  (v*k  fo  )  was  checked.  If,  in  fact, 

it  was  greater  than  JO  the  calculation  was  allowed  to  stand. 

However,  if  the  check  value  of  this  term  was  less  than  70, 

meaning  that  B  had  a  value  different  from  8.5,  then  the  new 

value  of  B  was  found  and  k  was  recalculated.  This  new  value  . 

s 

of  kg  was  then  accepted  without  further  checking. 

After  all  values  of  k  were  found  for  a  particular 

velocity  profile  (or  section)  the  results  were  inspected  and 

extraneous  values  were  rejected.  Often  such  values  were  those 

calculated  from  velocity  readings  near  the  surface,  or  near 

the  bed,  where  there  was  room  for  doubt  in  the  accuracy  of  the 

observed  reading.  The  remaining  values  of  k  for  that  profile 

s 

were  then  averaged.  The  average  values  of  all  three  sections 
for  each  depth  of  flow  were,  in  turn,  averaged  to  yield  a  final 
value  of  the  equivalent  sahdgrain  roughness  of  the  bed  form 
for  that  depth  of  flow. 

After  completing  these  calculations  for  all  depths 
of  flow,  the  results  were  converted  to  a  non-dimensional  form 
and  used  to  plot  the  desired  characteristic  curves  for  that 
particular  dune  shape. 

The  calculation  of  Chezy's  C  values  was  somewhat  more 
straight  forward.  For  each  depth  of  uniform  flow  the  recorded 
slope  and  discharge  were  used  in  equation  (2.2)  to  calculate 
the  required  value  of  Ghezy's  Q .  No  iterative  procedure  was  used. 
This  was  done  for  each  depth  of  flow,  and  each  dune  shape. 

Similarly,  the  values  of  C  were  then  converted  to 
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non  dimensional  form,  and  used  to  plot  the  required  curves  for 
that  dune  shape . 

The  observations  recorded  during  all  tests  are  pre¬ 
sented  in  Appendix  A.  These  readings  have  been  adjusted  for 
the  change  in  elevation  datum,  and  are  all  with  reference  to 
the  new  calculation  datum. 

Note  that  the  readings  recorded  on  the  laboratory 
data  sheets,  and  the  I.B.M.  punched  data  cards.,  on  file  in  the 
Department  of  Civil  Engineering,  are  not  adjusted,  and  are 
with  relation  to  the  observation  datum.  The  detailed  calcul¬ 
ation  sheets  for  each  velocity  profile  produced  by  the  computer 
also  on  file,  are  adjusted  to  the  new  datum,  however,  and  all 
readings  appearing  therein  are  with  reference  to  the  calcul¬ 
ation  datum. 

The  readings  shown  in  Appendix  A  a.re  all  identified 
as  to  test  run  by  showing  the  dune  shape  to  which  they  belong. 
Those  readings  taken  from  tests  run  on  sand  coated  dunes  are 
identified  by  an  R  after  the  dune  shape  in  the  table  heading. 
(For  example:  L/k  =  16  indicates  tests  done  on  the  smooth 
dune  shape  having  a  length  to  height  ratio  equal  to  16.  But 
L/k  =  16R  indicates  tests  done  on  the  same  dunes  after  extra 
roughness  in  the  form  of  a  sand  coating  had  been  added.) 


; 
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CHAPTER  IV 


EXPERIMENTAL  RESULTS 


4. 1  General 

This  chapter  contains  a  summary  of  the  experimental 
results  obtained  throughout  the  investigation.  The  actual 
data  collected,  and  the  results  calculated  from  that  data,  are 
tabulated  in  Appendix  A.  In  addition,  individual  plots  of 
these  results  for  each  test  run  are  contained  in  Appendix  B. 
The  actual  calculated  points  are  shown  on  these  graphs. 

The  graphs  contained  in  this  chapter  are  composite 
plots,  showing  all  related  tests  on  a  single'  sheet,  and  were 
obtained  by  merely  superimposing  the  appropriate  sheets  from 
Appendix  B.  Individual  calculated  points  are  not  shown,  only 
the  final  curves  are  reproduced. 

The  calculation  of  all  points  was  performed  as  out¬ 
lined  in  section  3.5.  The  graphs  contained  in  this  chapter, 
and  those  in  Appendix  B,  were  plotted  from  non-dimensional 
forms  of  these  results,  obtained  by  simple  arithmetic  division 
by  the  appropriate  factor  in  each  case.  These  calculations 
are  not  shown  here,  however  the  manner  of  conversion  from  the 
tabulated  results  to  the  plotted  non-dimensional  values  may 
be  easily  recognized.  In  most  cases  the  plotted  values  may 
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be  found  simply  by  inspection  of  the  calculated  results  in 

Appendix  A.  In  the  case  of  the  k  values.,  however,  each  plot- 

s 

ted  point  is  the  result  of  averaging  all  calculated  values, 
except  the  disregarded  ones,  for  a  given  depth. 

The  method  of  presentation  in  general  follows  that 

« 

outlined  in  section  2.4.  The  calculated  k  values  were  rendered 

s 

non-dimensional  by  dividing  them  by  k,  the  dune  height.  Thus, 
(ksA)  Sives  a  measure  of  the  effective  roughness  due  to  the 
dune  shape,  compared  to  the  actual  physical  height  of  the  dune. 
These  values  are  plotted  against  a  non-dimensional  depth,  d/k, 

found  in  a  similar  manner,  which  may  be  thought  of  as  an  indi- 

| 

cation  of  relative  roughness. 

Similarly,  the  Chezy's  C  values  were  rendered  non- 
dimensional  by  dividing  by  /gf.  Again,  this  may  be  thought  of 
as  a  measure  of  resistance  to  flow,  or  roughness.  These  values 
are  also  plotted  against  the  d/k,  relative  roughness,  parameter. 

4.2  Equivalent  Sand  Grain  Roughness 

The  non-dimensional  results  for  the  equivalent  sand 
grain  roughness  versus  depth  of  flow  have  been  shown  in  figs. 

4.1  and  4.2.  Fig.  4.1  deals  with  those  results  obtained  using 

smooth  concrete  bed  forms,  while  fig.  4.2  presents  those  results 

\ 

ootained  after  additional  roughness  had  been  added  in  the  form 
of  a  sand  coating. 


Each  of  these  graphs  contains  five  plots,  covering 
all  the  dune  shapes  tested.  The  dune  shape  for  each  curve  is 
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FIG.  4-2  VARIATION  IN  ROUGHNESS  WITH  FLOW  DEPTH  FOR  ROUGH  DUNES 
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indicated  on  the  graph.  The  purpose  of  both  these  figures  is 
to  indicate  the  dependence  of  roughness  on  flow  depth,  and, 
in  fact,  marked  trends  are  apparent. 

The  same  data  used  in  the  foregoing  plots  are  recon¬ 
structed  to  form  figs.  4.3  and  4.4.  Fig.  4.3  uses  the  data 
from  fig.  4.1,  but  it  is  plotted  in  a  different  manner.  Here 
the  roughness  kg//k  is  plotted  against  the  dune  shape  ratio, 
for  various  relative  depths  of  flow.  This  graph  is  intended 
to  point  out  the  dependence  of  roughness  on  the  shape  of  the 
bed  forms  which  are  present,  while  also  taking  account  of  the 
relative  depth  of  flow. 

Figure  4.3  deals  with  the  data  .obtained  in  the  tests 
run  on  the  smooth  concrete  dunes.  In  a  similar  manner,  fig. 
4.4  was  derived  from  fig.  4,2,  and  shows  the  results  of  tests 
run  on  the  same  bed  forms  after  the  additional  roughness  was 
added. 

The  effect  on  the  total  equivalent  sand  grain  rough¬ 
ness  of  the  addition  of  a  coating  of  uniform  sand  Is  shown  In 
fig.  4.5.  This  figure  was  derived  from  figs.  4.1  and  4.3,  and 
presents  the  Increase  in  roughness  caused  by  the  change  in  the 
surface  characteristics  of  the  dunes.  Two  ordinate  scales  are 

shown  on  this  graph.  The  first,  ^k  /k,  shows  the  change  in 

s 

the  kg/k  values  associated  with  each  relative  depth.  Plots 

are  included  for  each  dune  shape  tested.  The  second  ordinate 

scale  converts  the  change  in  k^/k  into  a  multiple  of  the  D 

s  50 

size  of  the  grains  used  in  the  sand  coating.  That  is,  it 
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demonstrates  the  effective  size  of  the  added  sand  grains  com¬ 
pared  with  their  actual  size. 

Figures  4.6  and  4.7  are  of  the  same  type  as  figs. 

4.1  and  4.2,  however,  they  were  arrived  at  in  a  somewhat 
different  way. 

Figures  4.1  and  4.2  were  plotted  from  values  calcul¬ 
ated,  as  was  previously  mentioned,  from  equation  (2.13). 

Thus  they  represent  values  of  roughness  derived  from  individ¬ 
ual  readings  of  velocity  at  given  elevations  in  the  flow. 
However,  figs.  4.6  and  4.7  were  derived  from  calculations 
based  on  equation  (2.14),  that  is,  on  the  observed  mean  veloc- 

I 

I 

ity  of  the  entire  flow. 

The  plots  of  figs.  4.6  and  4.7  ho  not  coincide  exactly 
with  their  counterparts  in  figs.  4.1  and  4.2.  This  should  be 
expected  due  to  their  different  methods  of  calculation,  and 
is  easily  explainable.  However,  a  similar  relationship 
between  roughness  and  relative  depth  is  apparent.  Again,  plots 
are  shown  for  each  dune  shape  tested. 

The  individual  plots  for  these  dune  shapes  are  not 
included  in  Appendix  B,  since  they  were  shown  only  for  purposes 
of  comparison.  Hence  the  original  points  calculated  are  shown 
on  these  graphs. 

Thus  the  relationship  between  the  non-dimensional 
equivalent  sand  grain  roughness, and  the  relative  depth  of  flow, 
and  the  dune  shape,  is  summarized  in  a  variety  of  forms  in 
figs.  4.1  through  4.7. 
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4 . 3  Chezy's  Roughness  Coefficient 

The  non-dimensional  values  of  Chezy1 s  C  were  treated 
in  exactly  the  same  manner  as  the  values  for  the  equivalent 
sand  grain  roughness. 

Figures  4.8  and  4.9  deal  with  the  results  for  the  . 
smooth  dunes ,  and  the  roughened  dunes  respectively.  The  non- 
dimensional  C//g"  term  is  a  measure  of  the  roughness,,  or  of 
the  resistance  to  flow.  It  is  shown  here  plotted  against  the 
non-dimensional  depth,  d/k,  which  could  he  called  a  relative 
depth,  or  a  measure  of  the  relative  roughness  of  the  boundary. 
These  plots  demonstrate  the  dependence  of  the  resistance  of 
the  bed  forms  on  the  depth  of  the  flow  in  the  channel.  Similar 
marked  trends  may  be  seen  in  figures  4.8  and  4.9.  In  each 
figure  curves  are  shown  for  all  dune  shapes  tested.  The 
Individual  calculated  points  are  not  shown  since  these  graphs 
are  meant  only  as  a  summary  of  results.  Individual  graphs  for 
each  dune  shape,  with  the  calculated  points  shown,  may  be  found 
In  Appendix  B. 

The  curves  shown  in  these  graphs  were  assumed  to  be 
straight  lines,  and  are  best-fit  lines  obtained  by  using  a 
least-squares  curve  fitting  technique. 

Figures  4.10  and  4.11,  dealing  with  results  from 
the  smooth  and  roughened  dunes  respectively,  show  the  depend- 
ance  of  the  resistance  to  flow  on  the  dune  shape,  for  various 
depths  of  flow.  These  figures  were  derived  from  figures  4.8 
and  4.9,  and  show  the  data  plotted  in  a  different  manner. 
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FIG„  4-8  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  FOR  SMOOTH  DUNES 
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F!G.  4-9  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  FOR  ROUGH  DUNES 
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FIG.  4-10  VARIATION  OF  ROUGHNESS  WITH  DUNE  SHAPE  FOR  SMOOTH  DUNES 
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VARIATION  OF  ROUGHNESS  WITH  DUNE  SHAPE  FOR  ROUGH  DUNES 


The  value  of  L/k  for  each  curve  is  indicated 
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In  figs.  4.10  and  4.11  the  non-dimensional  roughness,,  C/fg, 
is  plotted  against  the  height  to  length  ratio  of  the  dune 
profiles,,  k/L.  Curves  for  various  relative  depths  are  included. 

Figure  4.12  is  derived  from  figs.  4.8  and  4.9*  and 
shows,  the  change  in  the  resistance  to  flow  due  to  the  addition 
of  the  sand  coating  on  the  bed  forms.  Here  AC/fg,  the  change 
in  the  non-dimensional  resistance  coefficient  C//g^  is  plotted 
against  the  non-dimensional  depth  parameter  d/k.  Plots  are 
Included  for  each  dune  shape  tested.,  and  the  appropriate  length 
to  height  ratio  for  each  curve  is  indicated,  on  the  graph. 


CHAPTER  V 


DISCUSSION  OF  RESULTS 

5 • 1  General 

This  chapter  contains  a  general  discussion  of  the 
experimental  results  which  were  reported  in  Chapter  IV.  The 
form  of  the  resultant  graphs  is  discussed,,  and  possible 
sources  of  error  are  mentioned. 

5.2  Equivalent  Sand  Grain  Roughness 

Observation  of  figs.  4.1  and  4.2  shows  a  definite 
relationship  between  roughness  and  flow  depth.  It  is  seen 
that  an  increase  in  the  depth  of  flow  reduces  the  resistance 
due  to  the  bed  forms.  That  is,  a  smaller  relative  roughness 
produces  less  resistance  to  flow,  which  is  what  would  be 
expected. 

It  is  also  apparant  that  a  much  greater  change  is 
effected  in  the  case  of  bed  forms  having  shorter  wavelengths. 
However,  If  the  numerical  values  are  considered,  it  will  be 
noted  that  the  percentage  change  in  the  resistance  to  flow 
between,  for  example,  relative  depths  of  4.5  and  9.5,  though 
not  equal,  are  roughly  of  the  same  order. 

It  would  appear  that  the  resistance  of  the  dunes 
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with  shape  ratios  of  16  and  8  is  very  much  greater  than  that 
of  any  of  the  smoother  shapes.  The  reason  for  this  is  not 
immediately  apparent.  It  may  he  that  a  change  in  the  flow 
pattern  over  the  bed  forms  takes  place  between  the  curves  of 
the  dunes  having  shape  ratios  of  16  and  2k,  leading  to  larger 
.  eddies  and  resulting  in  greater  energy  dissipation. 

Examination  of  figs0  B1  to  BIO  in  Appendix  B  will 
show  that  the  scatter  of  points  from  which  these  curves  are 
drawn,  though  present,  is  not  by  any  means  exceptional.  It 
Is  felt  that  within  the  accuracy  of  the  experimental  tests, 

L 

and  for  the  original  purpose  of  the  study,  the  accuracy  of 
the  curves  is  satisfactory.  The  general  trends  of  the  results 
are  shown  clearly. 

It  will  be  seen  in  figs.  Bl,  B6,  b8,  and  B9  that 
the  point  or  points  at  the  shallow  end  of  the  curves  drop 
sharply  from  the  trend  Indicated  by  the  remainder  of  the 
points.  This  is  assumed  to  be  due  not  to  a  sudden  change  In 
flow  conditions,  but  to  the  fact  that  all  values  were  cal- 
culated  based  on  the  assumption  of  a  logarithmic  velocity 
profile.  It  is  assumed  that  at  these  small  depths  of  flow 
the  eddies  caused  by  the  flow  pattern  over  the  bed  forms  may 

disrupt  this  ideal  velocity  profile.  Yalin  (l4)  suggests  that 

V 

the  logarithmic  velocity  profile  may  not  be  valid  for  d/k 
values  of  less  than  6  or  7.  Thus,  these  values  were  assumed 
to  be  in  error,  and  were  not  used  in  drawing  the  final  curves. 

Figure  4.5.,  showing  the  change  in  roughness  caused 
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by  the  addition  of  the  sand  coating  to  the  bed  forms  is,  unfor¬ 
tunately,  rather  inconclusive.  In  general,  it  shows  that  the 
sand  coating  made  a  greater  difference  at  small  depths  of  flow, 
which  seems  reasonable.  However,  it  is  also  shown  that  the 
rougher  the  original  bed  form  pattern,  the  more  difference  that 
was  made  by  the  addition  of  sand  grains.  This  may  be  due  to 
the  greater  surface  area  of  the  rougher  patterns,  allowing 
more  flow  contact  with  the  sand  grain  roughness.  Again,  a 
rather  distinct  difference  is  noted  between  the  magnitude  of 
the  difference  found  for  the  16  and  8  shapes,  and  that  found 
for  the  flatter  shapes. 

The  fact  that  the  effective  increase  in  the  sand 
grain  roughness  is  greater  than  the  assumed  size  of  the 
coating  material  may  be  explained,  in  part,  by  the  unevenness 
of  the  sand  grain  layer.  Although  the  coating  was  made  reason¬ 
ably  uniform  some  irregularities  were  unavoidable.  Thus,  the 
effect  of  these  variations  was  magnified  beyond  the  size  of 
the  individual  sand  grains,  just  as  the  bed  forms  caused  an 
effective  roughness  greater  than  their  height. 

In  fig.  4.5  it  will  be  seen  that  all  plots  are 
roughly  parallel,  with  the  exception  of  that  for  dune  shape 
16.  This  deviation  is  unexplained,  and  is  likely  due  to  exper¬ 
imental  error.  The  same  accounts  for  the  sharp  upsweep  of  the 
L/k  =',32  curve  at  small  depths  of  flow. 

Although  the  effect  of  dune  shape  on  flow  resistance 
is  shown  in  the  previously  mentioned  figures,  it  is  more 
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apparent  in  figs .  4.3  and  4.4.  These  curves  demonstrate  that 
the  resistance  increases  sharply  as  the  wavelength  of  the  bed 
forms  is  decreased.  The  true  relation  between  roughness  and 
shape  may  not  be  well  represented,  however.  Unfortunately 
only  five  points  were  available  for  plotting  each  of  these 
curves.  In  addition.,  the  points  are  all  grouped  mostly 
toward  the  flatter  shapes.  Thus  a  small  change  in  the  values., 
for  instance,  of  L/k  =  16,  or  k/L  =  0.0625,  would  drasticaly 
change  the  shape  of  the  curves.  Since  the  point  derived  from 

the  L/k  =  8  tests,  or  k/L  =  0.125,  was  so  far  removed  from 

f 

the  remainder  of  the  points,  and  the  shape  of  the  curve  in 
this  interval  was  uncertain,  this  portion  of  the  curve  is  shown 
only  as  a  dashed  line.  The  general  trend  toward  increasing 
roughness  with  increasing  values  of  k/L  is,  however,  quite 
obvious.  The  reason  for  the  concave-up  shape  of  the  curves 
for  smooth  dunes  at  low  k/L,  while  the  corresponding  curves 
for  the  roughened  bed  forms  are  concave-down,  is  not  immed¬ 
iately  explainable.  Again,  this  is  possibly  due  to  insuffi¬ 
cient  data,  and/or  experimental  errors. 

Figures  4.6  and  4.7  present  the  roughness  calcul¬ 
ated  from  the  mean  velocity  readings.  It  will  be  seen  that 

theqe  curves  generally  give  higher  k  /k  values  than  figs.  4.1 

s 

and  4.2  respectively.  This  may  be  explained  by  considering 
that  figs.  4.1  and  4.2  were  derived  from  velocity  profiles 
taken  on  the  flume  centre  line,  largely  free  of  sidewall 
effects.  However,  the  figs.  4.6  and  4.7  are  based  on  the  flow 
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through  the  entire  cross  section,  which,  especially  in  a 
narrow  flume  such  as  the  one  used  in  this  study,  contains  a 
large  percentage  of  flow  which  is  retarded  due  to  boundary 
layer  growth  on  the  sidewalls  of  the  flume.  Thus  the  mean 
velocity  calculated  by  considering  the  entire  cross  section  • 
would  be  markedly  lower  than  that  observed  on  a  velocity 
profile  on  the  channel  centre  line.  Since  both  formulae 
(2.13)  and  (2.14)  are  similarly  derived,  and  make  no  allowance 
for  sidewall  effects,  the  roughness  or  resistance  calculated 
by  equation  (2.14)  is  bound  to  be  larger.  In  general,  both 
sets  of  curves  give  roughly  similar  values,  and  figs.  4.6 
and  4.7  provide  a  good  check  on  the  preceding  calculations. 

It  will  be  noted  that  in  both  figs  4.6  and  4.7  the  plot  for 
the  dune  shape  L/k  =  8  has  a  sharp  break  downward  at  small 
depths  of  flow.  This  is  again  assumed  to  be  due  to  the  lack 
of  a  logarithmic  velocity  distribution  at  these  small  depths 
over  a  very  rough  boundary. 

5.3  Chezy! s  Roughness  Coefficient 

The  results  of  the  Chezy's  C  calculations  are  in 
some  ways  more  satisfactory  than  those  of  the  equivalent 
sand  grain  roughness.  The  scatter  of  the  experimental  points 
was  less  severe,  as  may  be  seen  in  figs.  B-ll  to  B-20,  and  no 
untoward  breaks  were  seen  in  the  curves. 

Figures  4.8  and  4.9  demonstrate  the  effect  of  depth 
on  the  Chezy  parameter.  Again,  as  with  the  equivalent  sand 
grain  roughness,  an  unmistakable  trend  toward  decreased 
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resistance  to  flow  with  increased  depth  is  shown.  In  addition., 
it  would  seem  that  the  flatter  the  dune  shape,  the  more 
sensitive  it  is  to  changes  in  depth.  That  is,  there  is  a 
larger  change  In  roughness  for  the  dune  shape  having  L/k  =  40 
than  for  the  shape  L/k  =  8  when  the  depth  is  changed  from 
small  to  large  values. 

Again,  there  seems  to  be  a  distinct  break  between 
the  L/k  =  16  and  L/k  =  24  shapes,  similar  to  the  difference 
detected  in  the  ks  results.  Although  it  is  reasonable  to 
assume  that  if  the  difference  shows  up  in  one  roughness 
parameter  then  It  should  also  appear  in  the  other,  no 
satisfactory  explanation  can  be  made  for  this  phenomenon. 

In  fig.  4.8.  it  will  be  noted  that  the  results  for 
dune  shapes  having  ratios  of  40,  32  and  24  were  extended  to 
very  shallow  depths,  as  low  as  d/k  =  2.5.  These  were  among 
the  first  tests  performed,  and  the  practice  of  trying  such 
small  depths  was  later  discontinued.  This  action  was  taken 
due  to  the  extreme  difficulty  of  establishing  uniform  flow 
at  these  depths  with  the  equipment  used.  In  addition,  it  was 
decided  that  bed  form  heights  in  the  order  of  40/  of  the  depth 
flow  were  not  realistic  under  normal  circumstances. 

The  same  degree  of  confidence  should  not  be  placed 
in  these  readings,  taken  at  low  depths  of  flow,  as  those 
taken  in  deeper  water.  The  fact  that  the  L/k  =  24  and  the 
L/k  -  32  curves  cross  in  this  area  bears  out  these  suspicions. 

With  respect  to  accuracy  in  the  Investigation  it 
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should  be  pointed  out  that  In  fig.  4.8,  and  to  a  lesser  extent 
In  fig.  4.9,  the  general  appearance  of  the  set  of  curves  would 
be  Improved  if  the  L/k  =  24  curves  had  a  slightly  steeper 
slope.  Attempts  to  confirm  the  original  values  by  re-testing 
.  merely  duplicated  the  first  results.  However,  the  author  is 
not  convinced  that  the  curves  shown  are  correct,  and  it  is 
possible  that  experimental  error  caused  a  somewhat  flatter 
slope  than  should  actually  exist. 

The  difference  in  resistance  to  flow  caused  by  the 
addition  of  sand  coating  to  the  dunes  is  presented  in  fig.  4  12 

I  . 

As  with  fig.  4.5,  the  results  are  not  as  explicit  as  might  have 
been  hoped.  Examination  of  the  figure  will  show  that  the  change 
in  the  non-dimensional  Chezy  parameter,  C//g",  is  relatively 
insensitive  to  depth  of  flow.  That  is,  the  addition  of  the 
sand  coating  tended  to  add  a  certain  value  to  the  original 
roughness  coefficient,  and  this  added  amount  was  virtually  con¬ 
stant  with  changes  in  depth  of  flow. 

The  additional  roughness  was  not,  however,  constant 
over  varying  dune  shapes.  The  smoother  dunes  showed  a  much 
greater  change  in  resistance  to  flow  than  the  more  rugged 
profiles . 

The  fact  that  the  plot  for  L/k  =  32  is  not  parallel 
to  the  remainder  of  the  curves  is  not  immediately  explainable. 

It  is  probably  due  to  errors  in  the  investigation. 

The  relation  between  dune  shape  and  the  non-dimen¬ 
sional  Chezy  parameter  is  shorn  in  figs.  4,10  and  4.11.  As 
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with  figs.  4.3  and  4.4,  five  points  do  not  adequately  define 
the  underlying  relation.  Therefore  the  curves,  as  presented 
here,  may  not  be  correct  in  detail.  However,  it  would  seem 
that  an  unmistakable  trend  exists  toward  a  lower  C//g  value 
with  increased  k/L  values.  That  is,  more  resistance  to  flow  * 
with  a  more  rugged  bed  form  profile. 

It  is  apparent  that  if  the  G/fg  values  for  a  k/L  of 
0.04l6  ,  or  L/k  =  24,  could  be  lowered  slightly,  and/or  the 
C//g  values  for  a  k/L  of  0.0625,  or  L/k  =  16,  could  be  raised 
slightly,  then  a  smooth  concave-up  curve  would  result.  However, 
more  extensive  testing,  producing  more  points  to  plot,  would 

be  necessary  to  establish  whether  or  not  this  is  the  true 
relation. 

5.4  Possible  Sources  of  Error 

As  was  previously  mentioned,  one  possible  source  of 
error  in  the  experimental  procedure  was  the  application  of  a 
roughness  formula  based  on  a  logarithmic  velocity  profile.  A 
number  o±  the  velocity  profiles  taken  were  plotted  on  logarith¬ 
mic  paper  for  checking  purposes,  and  they  all  approximated  a 
straight  line.  However,  slight  deviations  did  occur,  especi- 
ally  near  the  bed  and  near  the  water  surface.  Most  of  the 
results  derived  from  doubtful  readings  were  rejected,  as  was 
mentioned  in  section  3.3;  however,  any  dubious  value  which 
was-  overlooked  would  lead  to  a  greater  scatter  in  the  final 
experimental  results. 
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Another  source  o±  error  which  may  be  considered,  is 
the  equation  (2.13).  This  equation  is  based  on  the  assumption 
of  fully  developed  turbulent  flow.  However,  during  the  cal¬ 
culation  of  some  of  the  results  it  was  found  that  the  parameter 
(v*ks/2v)  was  less  than  JO.  This .  indicates  that  fully  devel¬ 
oped  turbulent  flow  had  not  been  reached,  and  that  the  flow 
is  only  in  a  transitional  state.  In  a  few  cases  the  value 
of  this  term  dropped  below  5,  indicating  smooth  flow.  This 
condition  occurred  only  in  the  larger  d/k  region  of  the  smooth 
dunes  having  L/k  ratios  of  32  and  40. 

Under  these  conditions  the  viscous  sub-layer  theor- 

i 

etically  buries  all  roughness  elements.  That  is,  S  ,  the  thick¬ 
ness  of  the  viscous  sub-layer  is  greater  than  k  ,  the  equiv- 

s 

alent  sand  grain  roughness  height  of  the  roughness  elements. 

Thus  the  xlow  does  not  depend  on  the  roughness  characteristics 
of  the  boundary,  and  acts  as  if  the  boundary  were  smooth. 
Therefore  the  argument  could  be  made  that  equation- (2. 13) 
is  not  valid  for  these  calculations.  However,  rough  calcul¬ 
ations  would  show. thao  <$"  ,  at  its  thickest  value,  is  less  than 
one-tenth  the  actual  height  of  the  bed  forms.  Therefore, 

while  £  >  kg,  £<<k.  Thus  it  was  deemed  reasonable  to  proceed 

with  the  original  method  of  calculation. 

\ 

In  these  experiments  no  allowance  was  made  for  the 
resistance  to  flow  due  to  the  walls  of  the  flume.  All  resist¬ 
ance  was  attributed  to  the  bed  forms .  However,  since  the 
flume  walls  were  smooth  plexiglass,  and  contributed  much  less 


' 
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roughness  than  the  bed  forms.*  it  was  felt  that  the  resistance 
from  this  source  could  be  neglected. 

The  measurement  of  depths  which  determined  whether 
or  not  the  flow  was,  in  fact,  uniform,  was  another  possible 
source  of  error.  At  larger  depths  of  flow  the  water  surface’ 
was  glassy  and  the  depth  could  be  read  fairly  inccurately. 
however ,  when  ohe  flow  depth  was  reduced  to  2.5  inches  or 
less,  small  ripples  tended  to  appear  on  the  surface.  These 
tended  to  make  accurate  depth  determination  difficult,  so 
that  readings  could  be  taken  only  to  approximately  the  nearest 
one-tenth  of  an  inch.  Another  phenomenon  which  was  observed 
at  small  depths  oi  flow  was  a  dip  in  the  water  surface  over 
the  crest  of  each  dune.  This  is  normal  in  sub-critical  flow, 
but  is  not  usually  apparent.  The  small  depths  used  in  this 
study,  however,  magnified  this  effect.  In  order  to  avoid 
error  due  to  this  cause  all  depths  were  measured  at  similar 
points  on  the  bed  forms  at  each  end  of  the  test  section. 

Minor  errors  may  have  been  introduced  by  the  equip¬ 
ment  used  to  measure  the  discharge,  or  by  the  mini-flowmeter 
used  for  the  velocity  profile  determination.  However,  since 
both  these  systems  are  normally  quite  reliable  and  produce 
accurate  results  little  error  was  expected  from  these  sources. 

By  far  the  largest  errors  undoubtedly  arose  from 
measurements  of  slope.  The  scale  from  which  the  slope  was 
read  could  only  be  interpolated  accurately  to  the  nearest 
0.0002.  Any  closer  estimations  could  have  been  in  error. 
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POSSIBLE  PERCENTAGE  ERRORS  IN 
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Example:.  If'  error  in  the  observation  of  S  =  +10^ 

error  in  the  observation  of  y  =  -10% 
Then  the  resultant  error  in  k  =  +QA% 


TABLE  5 . 2 

POSSIBLE  PERCENTAGE  ERRORS  IN  C //g~ 
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Example:  If  the  error  in  the  observation  of  Q  =  -10% 

the  error  in  the  observation  of  d  =  +10% 
Then  the  resultant  error  in  C//g~  =-20. 2% 
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Since  in  many  cases  the  recorded  slope  was  of  this  order  of 
magnitude,,  it  is  possible  that  some  readings  could  have  been 
in  error  by  as  much  as  100$.  Any  such  errors  in  slope  esti¬ 
mation  would,  of  course,  effect  changes  in  the  calculated 
values  of  kg  and  C. 

During  all  tests  the  viscosity  of  the  water  was 
assumed  to  be  constant  at  the  value  for  70  F.  However,  the 
room  temperature  in  the  laboratory  probably  ranged  from  65  F. 
to  85  F.  during  the  course  of  the  study,  thereby  introducing 
a  slight  error.  However,  considering  the  magnitude  of  the 

errors  introduced  from  other  sources,  this  error  was  relatively 
unimportant . 

Tables  5*1  and  5*2  summarize  the  resulting  errors 
in  kg  and  C//g  respectively,  due  to  errors  in  the  observation 
of  various  parameters. 

5.5  Comparison  with  Results  of  Other  Studies 

A  number  of  studies  have  been  reported  in  the  lit¬ 
erature  which  have  used  mobile  bed  flumes  to  study  the  varia¬ 
tion  01  resistance  to  flow  with  changes  in  various  parameters. 
The  data  presented  in  these  references  may  be  compared  with 
the  results  obtained  in  this  investigation. 

Table  5 • 3  presents  data  already  published  in  various 
references.  In  each  case  the  original  data  were  analysed  to 
yield  C/C  gf  values  for  the  various  test  runs  performed,  based 
on  the  reported  values  of  slope,  depth  of  flow,  channel  width 


. 
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TABLE  5 . 3 

COMPARISON  WITH  OTHER  EXPERIMENTAL  RESULTS 


REFERENCE 

EXPERIMENTAL  RESULTS 

REF.  NO. 

DATA.  NO. 

c//F 

c/yr 

DIFF. 

5 

2-9 

8.14 

6.00 

+35.6 

9-1 

7.4o 

5.60 

+32.2 

9-2 

7.16 

5.60 

+27.8 

i 

9-3 

7.H 

5.70 

+24.7 

.» 

1 

9-4 

9.84 

6.00 

+64.0 

i 

i 

f 

9-5 

8.00 

5.50 

+45.5 

i 

9-6 

8.00 

6.00 

+33.3 

9-8 

6.90 

7.20 

-4.2 

9-9 

7.85 

5.75 

+36.6 

9-11 

9.20 

5.4o 

+70.4 

9-12 

9.65 

5.70 

+69.2 

9-15 

9.85 

6.10 

'+61.5 

6 

I 

12.00 

11.30 

+6.2 

II 

10.60 

11.30 

-6.2 

III 

9.60 

13.00 

-26.2 

IV 

10.80 

11.50 

-6.1 

7 

33 

11.90 

11.50 

+3.5 

la 

11.80 

11.40 

+3.5 

14 

12.30 

13.10 

-6.  r 

21 

12.60 

11.80 

+6 . 8 

19 

11.80 

11.60 

+1.7 . 

-e 

• 

TABLE  5*3  ( Coni’ d) 


72 


REFERENCE 

EXPERIMENTAL  RESULTS 

REF.  NO. 

DATA  NO. 

c//F 

c//g~ 

%  DIFF. 

7 

16a 

11.10 

10.80 

+2.8 

17a 

12'.  30 

11.90 

+3.4 

3a 

12.20 

11.60 

+5.2 

15 

2-8 

13.70 

13.80 

-0.7 

i 

2-11 

12.60 

11.20 

+16.2 

j 

i 

2-6 

10.40 

11.80 

-11.9 

f 

2-35 

11.50 

13.30 

-13.5 

i 

j 

2-9 

12.00 

11.10 

+8.1 

3-57 

9.50 

10.10 

-5.9 

3-54 

11.50 

14.00 

-17.8 

3-56 

10.60 

12.00 

-11.7 

3-43 

9.40 

11.00 

-14.5 

3-44 

11.30 

11.00 

*  +2.7 

4-33 

10.60 

11.20 

-5.4 

4-1 

10.70 

11.60 

-7.8 

4-14 

11.4o 

14.50 

-21.4 

4-2 

10.50 

10.80 

-2.8 

\ 

4-21 

11.20 

11.80 

-4.2 

4-19 

11.40 

11.70 

-2.6 

4-17 

11.40 

10.90 

+4.6 

4-3 

11.00 

11.70 

-6.0  ' 

4-34 

10.70 

15.90 

-32.7 

5-3 

11.20 

U 

10.70 

+4.7 

' 

. 
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TABLE  5*3  (Cont'd) 


REFERENCE 

EXPERIMENTAL  RESULTS 

REF.  NO. 

DATA  NO. 

C//1T 

c//T 

%  DIFF. 

15 

5-9 

10.50 

9.50 

+10.5 

5-1 

11.80 

12.10 

-2.5 

5-11 

9. 40 

7.00 

+34.1 

5-8 

9.40 

8.90 

+5.6 

5-6 

9. 40 

6.80 

-38.3 

5-12 

8.50 

6.90 

+23.1 

5-4 

12.80 

12.80 

0 

5-29 

11.80 

11.60 

+4.4 

6-l6 

10.50 

13.10 

-20.0 

6-17 

10.00 

10.20 

-2.0 

6-10 

11.20 

12.20 

-8.2 

6-7 

8. 40 

II.60 

-27.5 

6-8 

8.80 

12.60 

-30.0 

10-47 

15.90 

13.20 

+20.7 

10-43 

12.10 

13.00 

-6.9 

10-51 

7.50 

11.50 

-34.6 

21 

9.14 

_ 2^9_5 

-8.4 

' 
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and  mean  velocity.  Combining  these  values  with  equation  (2.2) 
yields  a  C//g  value  for  each  particular  run.  Then,,  using  the 
values  for  length  of  bed  form,  and  height  of  bed  form  also 
reported,  L/k,  and  d/k  ratios  were  calculated,  and  a  value  of 
C//g  was  obtained  from  the  experimental  results  of  this  study, 
as  shown  in  Chapter  IV. 

In  this  procedure  the  curves  of  fig.  4.8  were  used 
exclusively,  even  though  these  results  were  determined  for 
smooth  dunes.  By  using  the  data  reported  in  the  references 
the  D^o  size  of  the  sand  used  for  the  mobile  beds  was  found. 

In  addition,  the  average  length  of  dune  was  determined.  Com¬ 
bining  these  two  factors  gives  a  measure  of  the  scale  of  the 


surface  roughness,  compared  to  the  bed  forms.  The  average 
E^C)/L  value  each  study  was  as  follows: 

Reference  Average*  D^q/L 


6  0.000273 

7  0.000173 

15  0.000268 


21 


0.000485 


The  average  D^q/L  value  for.  this  present  Investiga¬ 
tion  was  0.00235.*  or5  about  five  to  fourteen  times  as  large 
as  the  values  for  the  reference  data.  In  view  of  these  results, 
it  was  assumed  that  the  dunes  referred  to  in  the  references 
could  be  considered  "smooth"  with  respect  to  the  sand  coated 
dunes  used  in  this  study,  and  that  a  better  comparison  could 
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TABLE  5.^-  (Cont’d) 
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be  made  with  the  results  of  the  tests  on  the  smooth  concrete 
bed  forms. 

The  comparison  of  the  experimental  results  with 
the  results  of  Vanoni  and  Hwang  (5)  are  not  as  close  as  those 
of  the  other  references.  In  their  study  Vanoni  and  Hwang 
arrived  at  a  dune  roughness  by  subtracting  a  grain  roughness, 
or  surface  roughness,  from  the  total  resistance  to  flow. 

This  surface  rougnness  was  determined  by  using  a  plot  of  the 
data  of  Nikuradse  for  sand  grain  roughness  in  pipes.  This  data 
was  obtained  for  smooth  walled  (i.e.  no  corrugations)  pipes, 
while  Vanoni  and  Hwang  had  a  surface  roughness  superimposed 
on  an  undulating  bed,  which  may  possibly  introduce  some  error 
and  account  for  the  larger  difference  between  their  results 
and  those  of  the  present  study. 

In  the  case  of  the  above  mentioned  study,  the  results 

were  compared  with  those  of  fig.  4.8  not  because  the  D^/L 

50' 

ratio  was  appropriately  small,  but  because  the  f"  factor 
calculated,  by  the  authors  should  show  only  the  roughness 
effects  due  to  the  form  drag. 

Table  5.4  presents  a  comparison  of  the  experimental 
results  with  values  calculated  by  proposed  formulae  by  which 
resistance  to  flow. may  be  predicted.  Formulae  put  forward  by 

Richa.rdson  (8),  and  Yalin  (.9)  are  examined.  The  formulae  are 

\ 

compared  with  results  from  both  smooth  and  rough  dunes,  of  which 
the  former  give  closer  agreement.  This  is,  as  before,  due  to 
the  scale  of  surface  roughness  in  these  tests,  compared  to  that 
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used  to  derive  the  formulae. 


5*8  The  Logarithmic  Velocity  Formula 

In  Chapter  II  the  following  expression  was  developed 
for  the  velocity  distribution  in  a  wide  open  channel: 


v  _  A  In  y  +  B  (2.1l) 

v*  ks 

By  transforming  natural  logarithms  into  logarithms 
to  the  base  ten.,  we  may  write: 


v  =  5-75  log  y  +  B 

v 


* 


k 


s 


Assuming  that  k^  and  B  are  constant,  we  can  say: 


v  _  A1  log  y  +  B* 
v*  k 

s 


(5.1) 


where:  A’  =  coefficient (usually  5.75) 

B'  =  constant 

The  relation  (5.1)  was  investigated  by  plotting 
values  of  v  against  log  y,  and  finding  the  best  straight  line 
fit  for  the  points.  The  value  of  the  coefficient  A’  was  then 
found  from  the  slope  of  the  curve. 

Unfortunately,  the  results  of  this  analysis  were 
not  conclusive.  A  summary  of  the  results  is  presented  below. 


and  the  error  with  respect  to  the  accepted  theoretical  value 
of  5.75  is  shown: 


Dune  Type 

A1 

average 

Error 

Smooth 

5.99 

+4.18^ 

Rough 

5.39 

-6.25^ 

Average 

5.69 

-l.o4 % 

■ 


b 


79 


A  number  of  investigators  (l),  (3),  (4),  (10),  have 
shown  that  some  variation  may  be  expected  in  this  coefficient. 
But  considerable  scatter  was  evident  in  the  results  of  the 
analysis  conducted  on  the  present  results ,  and  it  was  diffi¬ 
cult  to  define  any  firm  trends  in  the  results.  Individual 
values  of  A*  ranged  from  3.009  to  8.482,  indicating  a  varia¬ 
tion  of  Von  Barman's  coefficient  from  O.768  to  0.271.  This 
is  within  the  range  of  variation  discovered  by  Garde  and 
Paintal  (10)  for  alluvial  channels. 

It  was  originally  thought  that  a  slight  raising  or 
lowering  of  the  datum  might  lead  to  an  Improvement  In  the 
results  of  the  determination  of  A' .  However,  there  did  not 
seem  to  be  a  narrowing  of  the  scatter  of  A'  values  when  this 
was  attempted. 

The  complete  calculations  for  the  values  of  A'  for 
each  depth  and  dune  shape  are  not  shown  here,  but  are  Included 
with  the  remainder  of  the  thesis  data,  on  file  with  the 
Department  of  Civil  Engineering. 

It  will  be  seen  that  the. average  of  the  extreme 

values  of  A',  as  reported  above,  is  5.776,  or  close  to  the 

theoretical  value  of  5-75.  This  may  indicate  that  most  of 

the  variation  is  due  to  experimental  errors.  The  scatter  of 

the  values  calculated,  along  with  this  possibility,  showed 

that  the  assumption  of  a  value  of  A'  other  than  5.75  for 

purposes  of  calculating  k  was  unwarranted. 

s 

The  effect  of  changing  the  elevation  datum  was  studied 


80 


in  some  detail  by  calculating  values  of  A1  for  each  velocity 
profile  taken,  while  varying  the  elevation  datum  and  thus  all 
depth  measurements.  When  the  results,  for  each  dune  shape 
were  averaged  and  then  plotted  against  change  in  datum,  there 
was  no  improvement  in  the  scatter'  of  the  points,  and  it  was 
impossible  to  definitely  define  any  trend  in  the  results. 

For  this  reason  the  elevation  datum  was  left  at  its  original 
position,  midway  between  the  high  and  low  points  of  the  bed 
form  profile. 

1 

5.7  Comparison  with  Regime  Methods 

An  attempt  was  made  by,  Yalin  (16)  to  relate  the 
parameters  of  flow  and  particle  size  of  the  bed  material  to 
the  size  and  shape  of  the  bed  forms  which  may  be  expected  in 
a  channel.  The  relations  given  for  the  dimensions  of  the 
dunes  are  as  follows: 

k  _  3.  / 1  -  dcr  \  «  (5.2) 

d  5  \  —  / 

L  =  5d  (5.3) 

where :  d  =  depth  of  flow 

dcr  =  depth  corresponding  to 

competent  shear,  /X=J'S  d 
k  =  height  of  bed  form 
L  =  length  of  bed  form 

\  By  combining  these  equations  with  the  results  given 

in  fig.  4.9  it  should  then  be  possible  to  predict  the  channel 
resistance  if  the  flow  parameters  and  the  bed  material  size 
are  known. 

To  study  the  validity  of  this  approach  the  results 


.5 


' 

iVt  . 


8i 


are  compared  to  those  found  by  using  regime  methods. 

Table  5.5  shows  channel  geometry  and  Chezy's  C 
values  calculated.  1  or  various  discharges  by  the  regime  methods 
assuming  a  bed  material  diameter  of  0.718  mm., (the  same  size 
material  as  was  used  to  produce  the  roughness  coating  on  the’ 
smooth  bed  forms  in  this  investigation)  and  a  side  factor  of 
0 . 2_,  using  formulae  provided  by  Blench  (17).  In  addition, 
another  Chezy'  s  C  value,  to  be  called  C  ,  and  slope  have  been 

<u 

calculated  by  Yalin's  methods  and  are  Included  for  comparison. 
The  method  of  calculating  these  values  will  be  explained. 

Using  a  Shields^  diagram  for  competent  tractive  force, 
for  the  particle  size  used  in  this  study,  we  may  say: 

Scr  dcr  =  0.000133  (5.4) 

Using  equations  (5.2)  and  (5.4)  an  iterative 
solution  for  any  given  channel  may  be  found.  The  procedure 
used  is  as  follows: 


(1) 

(2) 

(3) 


Using  (5.2)  calculate  a  value  of  k. 


Using  d,  k,  equation  (5.3),  and  fig.  4.9, 
calculate  a  value  of  C  . 

u 


Substitute  Q,  b,  d 
to  calculate  S 
new  and  old  values 
the  procedure. 


,  and  C  into  Chezy's  equation 
If  the  difference  between  the 
of  S  is  significant,  continue 

u 


(4)  From  equation  (5.4)  calculate  a  new  d  and  hence 
a  new  value  of  k  for  step  (l).  cr 


Observation  of  Table  5.5  shows  that  the  agreement 
between  the  two  methods  is  not  very  good.  The  resistance  and 
slope  predicted  by  the  above  method  increases  with  increasing 
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discharge,  which  would  not  normally  be  expected. 

TABLE  5 . 5 

COMPARISON  OF  RESULTS  WITH  REGIME  METHOD 


Regime 

Method 

1 

i 

Yalin*  s 

Method 

Q 

d 

b 

S 

V 

c 

cy 

sy 

500 

3.4' 

63.4' 

.00022 

2.3  fps 

90 . 8  ! 

68.7 

.00039 

2000 

5.4 

126.8 

.00018 

2.9 

100.7' 

65.2 

.00042 

5000 

7.3 

200.4 

.00015 

3.4 

108.0 

64.1 

.00043 

Since  the  agreement  of  the  results  of  this  study 
with  those  derived  from  flume  experiments  is  reasonably  good., 
It  is  probable  that  equations  (5.2)  and  (5.3)  are  at  least 
partly  in  error. 

In  a  later  paper,  Yalin  (l8)  proposes  that  both 
sand  wave  height  and  length  can  be  expressed  as: 

/  2  V 

k,  L  =  f  /  Dv*  ,  v*  ,  d  ,  \ 

(nr  'jrjt  syj 

where:  =  unit  mass  of  sediment 

/o  =  unit  mass  of  fluid 

=  submerged  wTeight  of  sediment 

In  the  derivation  of  equations  (5.2)  and  (5.3)  the 
factor  d/b  was  not  considered.  This  omission  may  contribute 
some  of  the  error  between  the  two  sets  of  results  shown  in 
Table  5-5. 

The  term  (v*D/^)  ),  a  grain  size  Reynolds  number,  is 
mentioned  in  Yalin1 s  original  paper  as  a  criterion  for  the  type 
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O'  form  to  be  expected.  If  (v*D/^  )  <  20  ripples  will  form, 

but  if  (v^D/V;  )  ^  20  dunes  may  be  expected. 

White  and  MacMahon  ( 19 ) ^  and  Nordln  .(20)  also  crit¬ 
icize  Yalin's  results,  and  present  data  from  a  number  of 
sources  to  indicate  relationships  different  from  those  predicted 
by  equations  (5.2)  and  (5.3). 

Thus.,  while  Yalin’s  equations  give  results  accurate 
to  within  a  factor  of  two  or  three  when  compared  to  those 
predicted  by  the  regime  method.,  the  apparent  errors  may  be  due 
to  an  oversimplification  of  the  problem.  Research  on  the 
prediction  of  bed  profiles  is,  however,  as  yet  very  limited, 
and  hence  any  roughness  or  slope  values  predicted  by  the 
aforementioned  means  can  be  regarded  only  as  very  rough 
approximations . 
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CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  Conclusions 

A  review  of  the  current  literature  on  the  effect  of 
roughness  elements  on  open  channel  flow  showed  a  lack  of  knowl¬ 
edge  of  the  effects  of  shape  and  spacing  on  the  resistance  to 

■j 

flow.  This  was  particularly  evident  in  the  case  of  river 

I  , 

channel  bed  forms.  In  the  present  investigation,  the  use  of 
idealized  dune  shapes  in  a  laboratory  flume  presents  a  different 
approach  to  the  problem. 

The  tests  conducted  indicated  that  a  definite  rela¬ 
tionship  exists  between  resistance  to  flow  in  an  open  channel 
and  the  shape  and  height  of  the  bed  forms  which  occur  on  the 
channel  boundary.  In  addition,  the  grain  size  of  the  bed 
material,  or  the  surface  roughness  of  the  dunes,  was  shown 
to  have  a  considerable  effect  on  the  total  roughness. 

Examination  of  the  results  presented  suggests  a 
definite  increase  in  the  resistance  to  flow,  both  in  the  C//g" 
and  kg/k  parameters,  with  any  decrease  in  depth  for  a  given 
bed  form  pattern,  or  with  any  decrease  in  the  L/k  shape  ratio 
of  the  dunes  for  a  given  depth. 
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The  addition  of  a  surface  roughness  caused  increases 
in  the  total  equivalent  sand  grain  roughness  up  to  many  times 
the  size  of  the  actual  sand  grains  usedG 

The  results  obtained  in  this  investigation  were 
deemed  satisfactory  for  the  purposes  of  the  study.  The  methods 
used  proved  suitable  in  most  respects.,  and  ,  with  minor  changes, 
could  be  used  in  a  more  extensive  project. 

6.2  Recommendations 

It  is  recommended  that: 

(a)  a  more  extensive . series  of  tests  be  conducted 
utilizing  a  greater  range  of  dune  shapes,  and 
heights,  in  order  to  verify  or  limit  the  con¬ 
clusions  of  this  study. 

(b)  various  sizes  of  sand  be  employed  for  surface 
roughness  In  order  to  more  adequately  determine 
the  effect  of  this  factor  on  the  total  resist¬ 
ance  to  flow.  In  addition,  a  natural  river 
bed  gradation  of  sand  grains  should  be  tried, 
rather  than  a  uniform  size.  That  Is,  the  grain 
size  distribution  of  the  sand  used  should  follow 
a  log-normal  distribution. 

(c)  In  connection  with  (b)  above,  an  effort  should 

v  136  made  to  determine  the  best  characteristic 

size  of  the  bed  material  to  be  used  in  such 
plots  as  fig.  4.5,  showing  the  change  in  total 
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roughness  compared  to  the  actual  size  of  the 
bed  material. 

(d)  further  investigation  be  conducted  to  determine 
the  effect  of  the  breadth  to  depth  ratio  of  the 
flow  on  the  calculated  values  of  C/Vef,  and  k  /k. 
If  changes  occur  at  low  values  of  b/d,  then  it 
would  be  advisable  to  use  a  wider  flume  in 
future  studies.  This  would  permit  deeper  flows 
without  side  wall  effects,  and  would  allow 
extension  of  the  range  of  d/k  values. 

(e)  different  methods  of  measuring  the  depth  of  flow, 
and  slope  should  be  considered,  which  would 
allow  less  error  in  the  estimation  of  these 
values.  Greater  accuracy  in  slope  measurement 

is  especially  important. 

(f)  a  much  longer  flume  should  be  used  in  future 
studies  in  order  to  allow  accurate  confirmation 
of  the  presence  of  uniform  flow. 
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AiaTtKlSK  l;UT  USiQ  IiN  PLOTTING  TriF  CHARACTERISTIC  RPUGHNFSS  CURVES 


.V  A  K  K.  £  D  ,vl  h  AN  ASTCr.lSN  *fc*t  NOT  UStO  IN  PLOTTING  THE.  CHARACTERISTIC.  ROUGHNESS  CURVES 


ANALYSIS  OF  SAND  SAMPLE 


1)  Sieve  Analysis 

The  sand  used  to  form  the  additional  roughness 
coating  on  the  Led  forms  was  white  Ottawa  sand,  20-30  nominal 
size.  That  is,  the  grain  sizes  were  between  the  #20  and  #30 
sieve  sizes.  The  following  results  were  obtained  in  a  sieve 
analysis  of  a  sample  of  the  material: 

Material  larger  than  #20  2.4  gm.  0.126^ 

Material  smaller  than  #30  4.2  gm.  0.221^ 

Material  between  #20., #30  1899 •  8  gm  99.653^ 

Mesh  opening  in  #20  sieve  0.0331  inches 

Mesh  opening  in  #30  sieve  0.0232  inches 

Average  0.0281  inches 

2 )  Settlement  Analysis 

A  series  of  settlement  analyses  were  performed  on 
samples  of  the  roughness  coating  material,  using  a  Visual- 
Accumulation-Tube  Apparatus.  By  observation  of  the  recorder 
chart  produced  by  this  apparatus  the  grain  size  distribution 
of  a  sample  may  be  found. 

Some  difficulty  was  encountered  during  the  tests 
due  to  the  uniform  size  of  the  sample  grains.  It  was  observed 
that  the  sample,  when  released,  settled  virtually  as  a  unit. 
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making  it  impossible  to  accurately  plot  the  time-depth  curve 
of  the  sediment  accumulation.  The  results  obtained  were  as 
follows : 


Test 

1 

2 


0.570  mm. 


0.620  mm. 


Remarks 

Test  very  inaccurate 
-results  disregarded 


.3 


0.640  mm. 


Average 


O.63O  mm. 


3)  Grain  Measurement 

A  sample  of  the  sand  coating  material  was  placed 
under  a  microscope,,  and  the  grains  measured  against  the  eye¬ 
piece  grid.  Maximum  and  minimum  diameters  were  recorded  for 
twenty-five  grains,  and  the  average  diameter  calculated. 

Average  diameter  =  0.803  mm. 


4)  Average  Grain  Size 

To  summarize  the  results  of  both  the  tests: 

Test 

50 

Settlement.  O.63O  mm. 

Microscope  0.803  mm. 

Average  0.716  mm. 

Thus  the  average  diameter  according  to  these  tests 

was  O.716  mm.,  or  0.0282  inches. 

Referring  again  to  the  results  of  the  sieve  analysis, 
we  see  that  by  far  the  bulk  of  the  material  fell  between  the 
#20  and  #30  sieve  sizes,  which  have  a  mean  opening  of  0.0281 
inches.  This  compares  favorably  with  the  average  of  0.0282 
inches  determined  from  the  other  two  tests. 
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APPENDIX  B 


DETAILED 


EXPERIMENTAL 


RESULTS 
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FIG.  B-I  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k 
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FIG.  B-2  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  l6 
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FIG.  B-4  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH 
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FIG.  B-5  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  40 
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FIG.  B-6  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  8R 
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FIG.  B-7  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  l6R 
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FIG.  B-8  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  24R 
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FIG.  B-9  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  32R 
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FIG.  B-IO  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L./k  =  40R 
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FIG.  B-II  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k 
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FIG.  B-I2  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH 
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FIG.  B-I3  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k.  =  24 
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FIG.  B-I4  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  32 
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FIG0  B-I5  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =*  40 
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FIG.  B-I6  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  8R 
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FIG.  B-I7  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  l6R 
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-18  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  =  24R 
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FIG.  B-I9  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/ic 
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FIG,  B-20  VARIATION  OF  ROUGHNESS  WITH  FLOW  DEPTH  L/k  .=  40R 
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